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PREFACE 

This  work  represents  the  continuation  of  an  effort  begun  earlier 
with  the  preparation  of  the  Army  report  Fracture  Mechanics  Design  Hand- 
book. The  objectives  are  to  provide  an  Introduction  to  fracture  mechan- 
ics fundamentals  and  to  provide  a convenient  source  of  fracture  mechanics 
information  for  designers,  stress  analysts,  and  engineers.  The  first 
volume  was  limited  to  isotropic,  homogeneous  materials;  this  volume  is 
limited  to  filamentary  composite  materials. 

The  subject  of  fracture  mechanics  of  filamentary  composites  is  not 
yet  settled.  While  its  validity  for  making  strength  predictions  has 
been  confirmed  for  a number  of  composite  systems,  its  need  is  still  not 
clear.  The  need  is  not  clear  for  two  reasons.  While  it  is  true  that 
composites  are  subject  to  a certain  amount  of  crack-like  damage,  they 
are  not  generally  subject  to  fatigue  cracks  as  metals  are.  Also,  com- 
posites are  inherently  tough;  the  toughness  increases  with  the  unnotched 
strength.  The  question  of  need  will  perhaps  become  clearer  in  a few 
years  as  more  structures  are  built  of  composites.  The  experience  gained 
will  indicate  whether  brittle  fracture  and  the  need  to  control  it  are 
major  concerns. 

This  volume  contains  a summary  of  fracture  mechanics  techniques 
for  filamentary  composites.  An  elementary  introduction  to  the  notation 
and  constitutive  equations  of  composite  laminates  is  included  in  Chap- 
ter 1.  The  basic  stress  fields  on  which  fracture  mechanics  theory  is 
based  are  given  in  Chapter  2 along  with  several  numerical  solutions 
which  illustrate  certain  anisotropy  and  size  effects.  Chapter  3 is 
concerned  with  some  fundamental  ideas  of  fracture  behavior,  fracture 
mechanisms,  and  fracture  criteria.  Fracture  toughness  using  various 
specimen  geometries  for  epoxy  and  aluminum  matrix  composites  is  discussed 
in  Chapter  4.  Chapter  5 contains  two  theories  which  attempt  to  explain 
the  fracture  behavior  of  composites  with  stress  concentrations. 

For  a variety  of  reasons  a number  of  topics  related  to  the  present 
one  were  omitted  from  this  volume.  Chopped  fiber  composites  and  fatigue 
of  composites  were  not  included.  Because  work-of- fracture  measurements 
are  not  normally  specimen-independent,  work-of-fracture  ideas  were  not 
included.  They  are  primarily  used  to  obtain  the  relative  fracture  resis- 
tance of  different  materials.  Mode  111  type  crack  loading,  because  it 
it  is  believed  to  be  of  infrequent  importance,  was  not  mentioned. 

This  report  reflects  only  the  state-of-the-art  at  the  time  of  pub- 
lication and  it  is  expected  that  a stress  analyst  or  engineer  will  sup- 
plement the  data  presented  with  other  solutions  commonly  needed  in  his 
specific  work  as  they  become  available. 


No  pretense  of  originality  is  made  here.  Information  included  was 
gathered  from  many  sources  and  appreciation  is  extended  to  all  those  authors 
whose  work  is  cited. 


For  administrative  and  contractual  support  of  the  part  or  this 
project  performed  in  the  summer  of  1977,  the  authors  gratefully  acknow- 
ledge the  scientific  services  program  of  the  Army  Research  Office. 
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Chapter  1.  INTRODUCTION 


Fracture  mechanics  of  composites  is  a new  and  complex  subject. 
Reference  1 discussed  the  fundamentals  of  fracture  mechanics  of  iso- 
tropic materials;  the  following  chapters  discuss  the  fracture  mechanics 
of  composites;  its  analytical  foundation,  its  experimental  basis,  and 
the  fundamentals  of  its  application. 

1.1  Fracture  Mechanics  of  Composites 

The  concepts  of  linear  elastic  fracture  mechanics  were  developed 
primarily  for  applications  to  metals.  For  a number  of  reasons,  these 
concepts  are  currently  being  applied  to  composite  materials.  A strong 
advantage  of  filamentary  composites  is  that  they  can  be  tailor-made  with 
the  strength  and  stiffness  oriented  as  needed  for  a specific  application, 
so  that  high  strength-to-weight  ratios  result.  Because  of  this,  these 
composites  are  successfully  competing  with  metals  in  many  industries, 
especially  the  aerospace  industry.  In  material  comparisons,  fracture 
mechanics  provides  a consistent  way  of  ranking  the  fracture  resistance  of 
the  candidate  materials.  Furthermore,  the  fail-safe  design  philosophy 
of  damage  tolerance  or  damage  containment  requires  the  estimation  of 
residual  strength  of  damage  members.  This  has  encouraged  the  investi- 
gation of  fracture  mechanics  as  a tool  for  making  such  estimations. 

Fracture  mechanics  is  based  upon  the  occurrence  of  crack-like  flaws  in 
the  structure,  and  while  fatigue  cracks  normally  do  not  form  in  composites, 
flaw  damage  does  occur  in  a number  of  ways.  Ballistic  damage  to  military 
aircraft  is  an  outstanding  example.  Other  flaws  such  as  tool  scratches 
and  the  like  occur.  Some  flaws  are  inherent  to  the  production  process: 
internal  flaws  due  to  poor  wetting  of  the  fibers,  improper  adhesion  of 
the  matrix  and  fibers,  air  bubbles,  and  dirt  are  a few.  The  application 
of  fracture  mechanics  has  been  further  prompted  by  the  behavior  of  com- 
posite laminates  containing  notches  and  holes.  For  such  laminates, 
the  strength  does  not  conform  to  the  conventional  concept  of  a stress 
concentration  factor;  therefore,  fracture  mechanics  together  with  the 
idea  of  an  inherent  flaw  have  been  developed  for  making  strength  predic- 
tions in  such  cases. 

The  history  of  the  application  of  fracture  mechanics  of  composites 
is  quite  recent.  The  method  involves  determining  stresses  around  a 
flaw  which  is  modeled  as  a crack.  Although  heterogeneous  as  well  as 
anisotropic,  filamentary  composites  are  commonly  assumed  for  analysis 
purposes  to  be  representable  as  homogeneous  media.  Within  that  context 
Sih,  et  al.  [2,3]  provided  the  fundamental  solution  in  1965  for  the 
stresses  near  a flaw  tip  in  an  anisotropic  plate.  The  stresses,  although 
somewhat  different  from  the  isotropic  case,  retained  the  same  type  of 
stress  singularity  that  forms  the  basis  of  conventional  fracture  mechanics. 
Sih  defined  the  anisotropic  stress  intensity  factors  and  showed  for  certain 
fundamental  cases  that  these  were  the  same  as  their  isotropic 
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counterparts.  Sih  further  showed  that  these  stress  intensity  factors 
could  be  related  to  the  strain  energy  release  rate  of  Griffith  [4]. 

Wu's  experiments  [5,6]  then  showed  that  fracture  mechanics  concepts 
could  be  applied  to  composites.  His  experiments  were  limited  to  uni- 
directional composites,  with  the  initial  crack  parallel  to  the  fibers. 

In  a typical  composite  laminate  the  mode  of  fracture  is  far  more 
complex  than  it  was  for  Wu's  case.  For  angle  plies,  crack  extension 
is  frequently  accompanied  by  a general  area  of  crack  tip  damage:  delami- 
nation between  plies  and  splitting  between  fibers  within  plies.  Fracture 
surface  appearances  range  from  what  could  perhaps  be  described  as  smooth 
to  those  of  a "shaving  brush."  The  application  of  fracture  mechanics  to 
a wide  variety  of  laminate  configurations  been  demonstrated  by  Konish 
and  his  associates  [7,  8,  and  9]  through  fracture  toughness  testing 
programs . 

1.2  Introduction  to  Filamentary  Composites 

The  materials  under  discussion  here  are  limited  to  those  known 
as  filamentary  composites,  that  is,  materials  made  up  of  long  continuous 
fibers  embedded  in  a matrix  material  as  shown  in  Figure  1.  Common  reinforcing 
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Figure  1.  Fiber  reinforced  lamina. 
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fibers  are  glass,  graphite,  and  boron;  common  matrix  materials  are  epoxy 
and  aluminum.  The  fibers  have  extremely  small  diameters:  5.6  and  8 mils 

are  two  common  sizes  for  boron;  graphite  and  glass  are  an  order  of  magni- 

tude smaller,  having  diameters  of  0.3  to  0.4  mils.  The  amount  of  fiber 
material  is  normally  specified  by  giving  the  fiber  volume  ratio,  V/O, 
which  is  the  ratio  of  the  fiber  volume  to  the  total  volume. 

The  laminates  are  built  up  from  the  matrix  and  fiber  raw  materials 
by  a number  of  processes.  Shell  structures  are  typically  fabricated 
by  a filament  winding  process,  a method  whereby  a machine  winds  a con- 
tinuous bundle  of  the  fibers  wetted  with  the  epoxy  matrix  material  onto 

a mandrel  of  the  proper  shape.  Other  structural  laminates  are  fabricated 
by  the  lay-up  of  a tape  containing  both  the  fibers  and  epoxy  matrix 
materials.  The  3-in.  wide  tape  contains  the  untwisted  fibers  oriented 
in  the  longitudinal  direction.  The  successive  layers  of  tape  are  applied 
at  whatever  angle  is  required  to  give  the  resulting  laminate  its  desired 
properties  of  strength,  stiffness,  and  directionality.  A built-up  laminate 
is  shown  in  Figure  2.  After  the  laminate  is  built  up  from  the  tape  it 
is  subjected  to  a curing  process  of  heat  and  pressure. 
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Figure  2.  Filamentary  composite  laminate. 


In  the  literature  of  composites  a code  notation  is  used  to  desig- 
nate the  stacking  sequence  and  the  angular  direction  of  the  ply  fibers 
with  the  laminate  axis  (x-axis.  Figure  2).  A few  examples  [10]  will 
illustrate  the  method. 
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Laminate 


Code 


The  coding  sequence  starts  at  the  most  positive  z lamina.  The  "T" 
subscript  indicates  that  the  total  laminate  is  shown.  For  laminates 
having  mid-plane  symmetry,  a subscript  "sM  is  used  to  indicate  that  only 
half  the  laminate  is  shown,  as  follows: 

Laminate  Code 


[ 90/0^5}  b 


If  the  symmetric  laminate  has  an  odd  number  of  laminae,  the  code 
denoting  the  center  lamina  is  overlined,  indicating  that  half  of  the 
laminate  lies  on  either  side  of  that  lamina: 


i 


Laminate  Code 


In  informal  applications  of  the  code,  a laminate  may  be  referred 
to  by  giving  the  laminate's  set  but  not  the  number  of  sets.  For  example 
the  preceding  laminate  might  be  referred  to  as  simply  a (45/0/90) 
laminate.  8 

1.3  Constitutive  Equations 

The  constitutive  equations  are  thoroughly  discussed  [11];  the 
following  is  only  a brief  review.  First  an  anisotropic  sheet  in  a state 
of  plane  stress  is  considered.  The  in-plane  directions  are  denoted  as 
1 and  2.  Then  the  stress-strain  relationship  is 


where  the  stresses,  a^,  a and  and  the  strains  e^,  and  7j^»  are 

referred  to  the  1 and  2 axes.  The  quantities  are  the  elastic  con- 
stants comprising  the  stiffness  matrix.  The  stiffness  matrix  is  sym- 
metric so  that  six  independent  elastic  constants  completely  specify  the 
stress-strain  behavior. 


Consider  now  an  orthotropic  (three  mutually  perpendicular  planes 
of  symmetry)  sheet  as  represented  by  a lamina  of  filamentary  material. 
Figure  1.  It  is  assumed  that  the  lamina  is  homogeneous;  then  the  stress- 
strain  relation,  in  Equation  (1)  simplifies  to 
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where  the  elements  of  the  shiftness  matrix  are; 
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There  are  four  independent  elastic  constants;  E^  and  E^. 
Young's  moduli  In  the  1 and  2 directions,  respectively;  Gj^»  the  shear 
modulus;  and  the  major  Poisson's  ratio.  Note  also  that  from 

Equation  (3); 


V21E11  " V12E22 


The  compliance  matrix  can  be  found  by  inverting  the  stiffness 
matrix.  Results  are; 
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where  the  elements  of  the  compliance  matrix  are. 
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It  is  important  to  note  that  the  stress-strain  relations  in  Equations  (2) 
and  (4)  are  for  a specially  orthotropic  lamina,  where  the  constitutive 
relations  are  referred  to  the  principal  axes  (1,2)  of  the  orthotropic 
lamina. 


Normally  the  principal  axes  of  the  lamina  do  not  coincide  with  the 
reference  x and  y axes  for  the  laminate.  In  that  case  the  constitutive 
relationship  must  be  transformed  through  an  angle  9 to  the  laminate  axis 
x and  y as  shown  in  Figure  1.  For  brevity,  the  details  which  are 
included  in  Reference  11  will  be  skipped.  The  resulting  stress  strain 
relation  is; 
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where  Q^,  the  components  of  the  stiffness  matrix  referred  to  the 
arbitrary  x-y  axis,  are  [11]: 


Qjj  18  Qj^  cos4  0 + 2(Q12  + 2 Q66)  sin2  0 cos2  0 + Q22  sin4  0 
Q22  “ Qu  sin^  9 + 2(Q12  + ^ Q^6)  sin2  0 cos2  0 + Q22  cos4  0 

Ql2  = (Q1X  + Q22  " 4 8in2  9 co®2  9 + Q12  (8in4  6 + co®4  0> 

%6  = «ll  + ^22  - 2 Q12  " 2 <W  8in2  9 COs2  9 

+ Q66  (sin4  0 + cos4  0) 

Ql6  = (Q11  " Q12  - 2 Q66>  8in  9 COs3  9 + (Q12  ' Q22 

+ 2 Q66)  sin3  0 cos  0 

^26  * «11  - ^12  - 2 %6>  8l”3  9 c°"  9 + «12  - <J22 

3 

+ 2 Q, , ) sin  0 cos  0 

DO 


None  of  the  terms  in  the  stiffness  matrix  are  zero  now.  In  fact, 
the  stiffness  matrix  has  the  same  form  as  for  a fully  anisotropic  plate. 
Equation  (1).  However,  there  are  still  only  four  independent  elastic 
constants.  The  terms  7$^  and  are  merely  linear  combinations  of  the 

first  four  Q's.  Of  course,  this  stiffness^  matrix  can  be  inverted  to 
find  the  corresponding  compliance  matrix  S for  an  orthotropic  plate. 
Thus , 
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where  [S]  is  the  inverse  of  [Q],  or  [S] 


1.4  Laminate  Constitutive  Equations 

Given  the  location  and  orientation  of  all  the  plies,  it  is 
possible  to  compute  the  macroscopic  elastic  properties  of  a laminate 
from  the  elastic  properties  of  its  constituent  plies.  To  do  so,  the 
constitutive  equations  for  each  k ply  shown  in  Figure  3,  are  first 
referred  to  the  laminate  axes  x and  y,  then  the  stress  resultants  are 
computed  across  the  laminate  thickness,  B,  by  integrating  the  stresses 
across  each  lamina  thickness. 
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Figure  3.  Notation  for  lamina  coordinate  within  a laminate. 


In  general,  the  results  are  rather  lengthy:  the  in-plane  stress 
resultants  not  only  depend  upon  the  mid-plane  strains  but  upon  the 
bending  curvature  as  well;  the  laminate  exhibits  coupling  of  extension 
and  bending.  This  means  that  applied  extension  causes  bending  or 
warping  or,  conversely,  an  applied  moment  causes  in-plane  forces.  This 
is  normally  an  undesirable  situation  in  that  panels  subjected  to  tension 
would  warp  out  of  plane  even  if  the  tension  were  due  to  thermal  stresses 
resulting  from  a change  in  temperature. 


It  can  be  shown  that  the  bending-moment  coupling  disappears  if  the 
laminates  possess  mid-plane  symmetry  [ 11] . Mid-plane  symmetry  means 
that  for  every  ply  above  the  laminate  mid-plane  there  is  one  of  equal 
properties  and  orientation  an  equal  distance  below  the  mid- plane.  Most 
laminates  have  mid-plane  symmetry  to  prevent  the  undesirable  warping 
mentioned  previously.  Consequently,  all  of  the  following  will  be  restricted 
to  mid-plane  symmetric  laminates.  Moreover,  applied  moments  are  not 
considered;  therefore,  the  moment -curvature  constitutive  relations  will 
not  be  required.  With  these  limitations  then,  the  stress  resultants 


(N  , N , N ) on  the  laminate  element  in  Figure  3 are  related  to  the 
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where 


k?i  (QiJ>k  (hfc  " hk-l) 


(10) 


Thus  for  mid-plane  symmetry  the  in-plane  laminate  stiffness  is  essentially 
the  sum  of  the  stiffnesses  of  the  constituent  plies  Q.  , all  referred  to 
the  common  laminate  axes,  x and  y.  J 


If  one  works  in  terms  of  the  average  laminate  stresses,  a,  a , t , 
where  - xx  yy  xy 
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then  the  laminate  stress-strain  relationship  becomes 


cr 
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(12) 


where 


[C]  = ^ [A] 


(13) 


The  compliance  matrix  [a]  can  be  found  by  taking  the  inverse  of  the 
stiffness  matrix,  [ C] . That  is, 

[a]  - [Cl”1  (14) 

so  that 


(15) 


Equations  (12)  and  (15)  are  the  laminate  constitutive  equations.  The 
elements  C^,  of  the  stiffness  matrix  are  computed  from  Equations  (10) 
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and  (13);  and  Che  elements  a,  , 
Equation  (14).  ^ 


of  the  compliance  matrix  are  given  by 


If  for  every 
properties  with  a 


ply  having  a +9 
- 8 orientation. 


orientation  there  is  a ply  of  equal 
then  a^  » a ^ " 0 and  the  laminate 


is  specially  orthotropic,  having  only  four  independent  elastic  constants. 
In  that  case,  the  elements  a^  of  the  compliance  matrix  are  defined  in 

terms  of  the  laminate's  engineering  constants  in  a manner  identical  to 
Equation  (5) , which  were  derived  for  a single  orthotropic  lamina. 

A laminate  which  contains  an  equal  number  of  laminae  of  +9  and  -Q  fiber 
orientation  is  referred  to  as  a "balanced"  laminate. 


The  term  quasi-isotropic,  sometimes  used  in  the  literature,  refers 

to  a particular  type  of  laminate  which  is  nearly  isotropic  in  stiffness. 

Examples  are:  a three-ply  laminate  with  a [0°/±60]  orientation,  and 

8 

a four-ply  laminate  with  a [0°/±45°/90#]  orientation. 

8 
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Chapter  2.  CRACK  TIP  STRESSES  FOR  ANISTROPIC  MATERIALS 


l 


1 


The  application  of  fracture  mechanics  is  based  upon  the  mathematical 
form  of  the  stresses  around  the  tip  of  a sharp  flaw.  Results  leading 
to  basic  definitions  have  been  obtained  for  fundamental  geometries  in 
anisotropic  materials.  Some  of  these  results  will  be  discussed.  Addi- 
tionally, some  numerical  solutions  which  demonstrate  certain  geometric 
and  material  effects  for  selected  filamentary  laminates  will  be  consid- 
ered. 

2.1  Crack-Tip  Stress  Fields  in  Anisotropic  Plates 

In  Chapter  1,  the  form  of  the  stress-strain  relationship  for  an 
anisotropic  plate  (or  laminate)  was  noted  to  take  the  form 


re  ~\ 

XX 

”all  al2  a16 

‘‘xx' 

€yy 

a 

a12  a22  a26 

ayy 

7*y. 

a16  a26  a66_ 

.V. 

where  the  a^  are 
anisotropic  plate 


the  compliance  coefficients  for  the  plate.  If  the 

is  also  a laminated  plate  then  the  stresses  a , a , 

xx  yy' 


and  t must  be  considered  to  be  stresses  averaged  over  the  laminate 
xy 

thickness.  Here  the  bar  denoting  average  has  been  dropped  and  the 
stresses,  when  applied  to  a laminate,  will  be  understood  to  represent 
average  stresses. 


The  form  of  the  crack-tip  stresses  for  an  anisotropic  plate  shown 
in  Figure  4,  has  been  presented  by  Sih  et  al.  [2]  and  others  [3,12] 
using  either  those  methods  discussed  by  Lekhnitskii  [ 13]  or  the  Westergaard 
stress  function  approach  [ 12] . Use  of  the  equilibrium  equations  together 
with  the  compatibility  equations  results  in  a fourth-order  governing 
equation.  Solving  this  equation  involves  finding  the  roots  ^ of  the 
characteristic  equation, 

aH  ^ " 2 ai6  ^ + <2  a12  + a66^2  " 2 “26  + a22  " 0 * <17) 


These  roots  are  complex  or  purely  imaginary  and  occur  in  conjugate 
pairs,  i.e.,  and  ^ * where  the  bar  denotes  conjugate. 

Omitting  the  details,  the  stresses  and,  displacements  [2]  are  given  in 
the  following  equations  for  two  cases:  (1)  the  loads  symmetric  about 
the  crack  plane  and  (2)  the  loads  skew  symmetric  about  it. 


< 
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Loads  8kew-synmetric: 


a " — T Re 

“ J2m 

hi 

a  *  **  — — Re 

**  J2^ 

Kl1 

Txy  = 7=^  Re 
y Jim 


u»  sin  0 


Lt,  sin  0 


cos  0 + 4-  sin  9 Jco8  9 + 4,  sin  9 


cos  0+4,  sin  0 ycos  0 + 4-  sin  0 


u - %/?  Re 

V ^ I \“2y  cos  0+^2  sin  ^ ’iy  cos  0 + sin  0 

The  polar  coordinates,  r and  0,  shown  in  Figure  4;  ,4^  and  42  are 
roots  of  the  characteristic  equation,  and  are  stress  intensity 

factors,  u and  v are  displacements  in  the  x and  y direction,  respectively, 
and  the  quantities  p^  and  q^  are  convenient  combinations  of  certain  elastic 

properties  given  by: 


an  “j  + ai 


*16  **) 


*12  “j  +77  ' “26 


(J  - 1,2) 


From  Equations  (18)  and  (19),  it  can  be  noted  that  the  stresses 

-1/2 

exhibit  a crack-tip  stress  singularity  of  r ¥ as  in  the  isotropic  case. 
The  angular  distribution  of  the  stresses,  that  is,  the  variation  of  the 
stresses  with  0,  depends  upon  the  material  properties  through  the  roots 
4^.  The  displacements  also  exhibit  an  anisotropic  effect:  for  sym- 
metric loading,  displacements  are  not  necessarily  symmetric;  for  skew- 
symmetric  loading,  displacements  are  not  necessarily  skew- symmetric  as 
in  the  isotropic  case.  This  means  that  skew- symmetric  loading  causes 
crack  opening  in  addition  to  forward  sliding;  symmetric  loading  causes 
crack  forward  sliding  in  addition  to  crack  opening.  This  displacement 
effect  disappears  if  the  crack  is  oriented  along  one  of  the  material's 
principal  directions. 


The  stress  intensity  factors,  Kj.  and  K^,  represent  the  strength 

of  the  stress  singularities.  They  have  been  defined  in  a way  which  is 
analogous  to  the  isotropic  case.  Formulas  for  finding  and  from 

the  appropriate  stress  functions  are  given  in  Reference  2.  Omitting 
those  details,  consider  the  results  for  some  fundamental  cases.  For 
example,  the  results  for  the  geometry  shown  in  Figure  5,  a crack  of 
length  2a  in  an  infinite  plate,  are: 

1 4-0  («a)1/2 


Figure  5.  Tunnel  or  line  crack  in  an  anisotropic  plate 
loaded  with  biaxial  tension. 

Results  for  the  probelm  in  Figure  6 are; 


Kj!  " T (Jt«r 
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Figure  6.  Tunnel  or  line  crack  in  an  anisotropic  plate 
under  remote  shear. 

As  a final  example,  the  results  for  the  wedging  force  problem  in 
Figure  7 are: 

mrsT/i 

(jra) 


(23) 


•4i 


- o 


Figure  7.  Tunnel  or  line  crack  with  a crack  surface  wedging  force. 


It  will  be  noted  that  the  expressions  for  Rj.  and  R^,  Equations  (21), 

(22),  and  (23)  are  identical  with  their  isotropic  counterparts.  They 
depend  upon  the  loading  and  crack  length  the  same  way  as  for  the  iso- 
tropic case.  For  cracks  in  Infinite  anisotropic  bodies  where  the  stress 
resultant  on  the  crack  is  zero,  these  results  are  true  in  general:  the 
anisotropic  stress  intensity  factor  is  the  same  as  for  the  isotropic 
case  with  the  same  geometry  and  loading.  In  other  words,  the  stress 
Intensity  factors  are  not  material  dependent  for  cracks  in  infinite 
plates.  If  the  stress  resultant  on  the  crack  surface  is  not  zero  or  if 
stress  free  boundaries  occur  near  the  crack,  then  the  anisotropic  stress 
Intensity  factor  will  differ  from  the  isotropic  one.  ThiB  latter  case 
will  be  discussed  in  more  detail  later. 


One  must  be  careful  in  using  formulas  for  stress  intensity  factors. 

Two  similar  definitions  of  both  Kj  and  differing  by  the  factor 

sA,  are  in  common  usage.  The  two  definitions  are  normally  distinguished 

by  the  use  of  lower-case  or  upper-case  K;  this  is  the  notation  adopted 

here.  That  is,  Kj  «,/«  kj  and  ■v/7kII. 

The  basis  of  the  equivalence  of  the  Griffith  concept  of  fracture 
and  the  stress  intensity  approach  is  the  relationshup  between  the  stress 
intensity  factor  and  g,  the  strain  energy  release  rate.  As  in  the  iso- 
tropic case,  the  strain  energy  release  rate  can  be  determined  by  calcu- 
lating the  work  done  by  the  crack-tip  stresses  in  closing  the  crack  over 
a vanishingly  small  distance  at  the  crack  tip.  To  avoid  almost  insur- 
mountable mathematical  complications  in  performing  this  calculation,  it 
is  assumed  that  the  crack  extends  in  a collinear  manner  although  in  an 
anisotropic  material  this  frequently  will  not  be  the  case.  If  symmetric 
loading  (Mode  I)  and  skew-symmetric  loading  (Mode  II)  are  considered 
separately,  results  are  [2], 


a22 


all 


+ ^2 

. ^ . 
[4i  + u2] 


where  g^  and  g^  are  the  strain  energy  release  rates  for  symmetric  and 

skew- symmetric  loading,  respectively.  When  both  modes  are  present,  then 
the  cross-inf luences  of  the  displacements  cause  cross-product  terms  of 
and  in  the  expressions  for  g,  namely. 


*1  fc  S + ^2>  + 'll! 
gI  “ " T a22  H T- 

L V2  J 

*11  fci  *1  1 

«II  " T"  an  Hyf  (^1  + ^2>  +y= 


where 


*total 


" ®I  + *11 


For  the  special  case  of  an  orthotropic  material  with  the  crack  on  one 
of  the  material1 8 planes  of  symmetry,  the  two  modes  are  independent,  so 
that 


i 


In  summary,  it  is  well  to  note  the  numerous  similarities  between  the 
isotropic  and  anisotropic  cases.  Equations  (18)  and  (19)  indicate  that 
for  a given  material  the  crack-tip  stress  distribution  for  different 

specimen  geometries  will  be  identical,  and  the  r~1^2  stress  singularity  of 
isotropic  materials  is  preserved.  Furthermore,  the  stress  intensity  fac- 
tor depends  upon  the  load  and  specimen  geometry  in  the  same  way  as  for 
isotropic  problems,  and  is  related  to  Griffith's  strain  energy  release 
rate.  In  fact,  as  already  noted,  the  stress  intensity  factor  for  a 
crack  in  an  anisotropic  plate,  except  for  two  conditions,  is  exactly  the 
same  as  for  a crack  in  an  isotropic  plate  having  the  same  loading  and 
geometry.  The  two  conditions  are  that  the  plate  must  be  infinite  in 
extent  and  the  resultant  of  any  loads  on  the  crack  surfaces  must  be  zero. 
The  finite  plate  effect  is  a detail  which  will  be  dealt  with  in  the  next 
article.  It  thus  seems  reasonable  to  assume  that  the  stress  intensity 
factor  concept  can  be  employed  to  predict  fracture  for  anisotropic  mate- 
rials as  is  done  for  isotropic  materials.  Proof  that  this  idea  also 
applies  to  filamentary  composites,  which  are  heterogeneous  as  well  as 
anisotropic,  must  be  based  upon  experimental  investigations  of  crack 
extension  in  various  composites. 

2.2  Effect  of  Finite  Specimen  Size 

To  apply  conventional  fracture  mechanics  to  composite  materials  a 
knowledge  of  the  stress  intensity  factor  is  required.  This  is  needed 
both  for  cracked  structural  components  under  analysis  for  residual 
strength  and  for  specimens  proposed  for  fracture  toughness  testing.  A 
considerable  collection  of  stress  intensity  factor  solutions  already 
exists  for  isotropic  materials  of  various  geometries  and  loading  condi- 
tions [1,  14].  Thus,  it  would  be  advantageous  to  understand  to  what 
extent  isotropic  stress  intensity  factors  could  be  applied  to  anisotropic 
cases  to  recognize  problems  for  which  this  practice  would  be  expected  to 
be  accurate  and  to  have  a basis  for  judging  the  degree  of  accuracy.  As 
mentioned  before,  for  problems  involving  free  boundaries  near  the  crack 
tip,  the  stress  intensity  factor  exhibits  a dependency  upon  the 
anisotropic  material  properties.  Solutions  for  such  problems  involve 
numerical  techniques  such  as  the  boundary  integral  equation  method  [15, 

16,  17]  finite  elements  [18],  and  boundary  collocation  [19]. 
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The  double-edge-cracked  (DEC)  and  the  center- cracked  tension  (CCT) 
specimens  are  familiar  fracture  toughness  specimens  for  isotropic  mate- 
rials. Thus  it  is  not  suprising  to  find  that  numerical  work  for  com- 
posite laminates  has  been  largely  devoted  to  these  specimens. 

Snyder  and  Cruse  [ 16]  determined  the  finite  correction  factor  for 
the  CCT  specimen  for  a number  of  practical  laminate  configurations.  The 
lamina  properties  were  chosen  as 


ER  - 21  X 106  psi 
E22  = 1.7  X 106  psi 

G12  “ 1,4  X 1()6  p8i 


(27) 


which  are  representative  of  graphite-epoxy.  A sample  of  the  results 
are  shown  in  Table  1 and  Figure  8.  The  lamina  angles  are  relative  to 
the  crack  direction.  All  results  are  for  a plate  length  to  width  aspect 
ratio  of  2L/W  = 3.  Results  are  given  in  the  form. 


where  Y is  the  finite  plate  width  correction,  K^,  is  the  Mode  I stress 

intensity  factor,  o is  the  applied  stress,  and  a is  the  half  crack 

length.  The  value  of  Y for  a crack  in  a infinite  isotropic  plate  is  ^/jT. 

In  Table  1 the  results  for  the  Isotropic  case  should  first  be  compared 

with  the  quasi-isotropic  laminate,  (0/±45/90)  . Except  for  2a/W  = 0.8, 

s 

the  results  for  the  two  are  within  a fraction  of  17.  of  each  other; 
therefore, the  isotropic  and  quasi-istropic  results  agree  well.  This  must 
be  kept  in  mind  in  comparing  the  results  on  Figure  8,  Note  that  the 
amount  of  departure  of  the  angle-ply  results  from  those  for  the  quasi- 
isotropic  case  is  a measure  of  anisotropy.  All  results  are  fairly  close 
to  Isotropic  (quasi-isotropic)  results;  with  the  greater t departure 
occurring  for  the  (±45)  laminate.  The  results  for  the  (±45)  laminate 

8 B 

are  approximately  117.  higher  than  for  the  isotropic  plate.  The  stress 
intensity  for  the  (90)  laminate  (fibers  perpendicular  to  the  crack) 

is  nearly  the  same  as  for  the  isotropic  case;  the  results  for  the  (0)g 

laminate  (fibers  parallel  to  the  crack)  are  slightly  lower  than  for  the 

isotropic  case.  The  variation  of  the  stress  intensity  with  ply  angle, 

a,  for  the  (±a)  laminates  is  shown  in  Figure  9 for  2a/W  “ 0.6. 
s 
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TABLE  1.  FINITE  WIDTH  CORRECTION  FACTORS  FOR  CENTER-CRACKED 
PLATES  (2L/W  - 3)  f 161 


Lamina  angles  measured  from  crack 
axis. 

Constant  displacement  applied  to 
ends  of  specimen. 

See  Figure  8. 
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Again,  the  greatest  departure  from  the  isotropic  (quasi-isotropic)  case 
is  for  a - 45°.  Snyder  and  Cruse  [16]  used  both  a constant  displacement 
loading  and  a constant  stress  loading  on  the  ends  of  the  plate.  The  two 
cases  were  within  2%  except  for  the  (90)  laminate  where  the  constant 

displacement  loading  was  approximately  37«  lower  than  the  constant  stress 
loading.  For  a shorter  plate  this  difference  would  increase. 

Konish  [ 20]  using  the  boundary  integral  technique,  obtained  results 
similar  to  those  in  Reference  16  for  the  CCT.  His  results  are  written 
in  the  form 

ICj.  - HY*  Oy/a  (29) 


where  Y*  is  the  isotropic  geometry  correction  factor  which  is  a function 
of  the  specimen  geometry  only  and  H is  an  anistrophy  factor.  The  ratio 
of  the  anisotropic  stress  intensity  factor  and  its  isotropic  counterpart, 
H,  is  defined  by 


H 


*1 

y*otT 


K_:  anisotropic 

K^:  isotropic 


(30) 


Thus,  the  deviation  of  H from  unity  is  a direct  measure  of  the  effect  of 
material  anisotropy  on  the  stress  intensity  factor. 


The  variation  of  the  anisotropy  factor  H with  crack  length  and  ply 
angle  is  tabulated  in  Table  2 and  shown  graphically  in  Figure  10  for  the 
CCT  specimen.  Results  are  likewise  given  for  the  DEC  specimen  in  Table  3 
and  Figure  11.  For  both  specimens  the  properties  of  the  graphite  plies 
were 


En  - 20.5  X 106  psi 
E22  - 1.37  X 106  psi 
G12  - 0.752  X 106  psi 


(31) 


v 


12 


0.31 


> 


which,  except  for  G^,  are  very  nearly  the  same  as  those  used  by  Snyder 

and  Cruse  [ 16] , Equations  (27) . The  ply  angle  a in  Tables  2 and  3 and 
in  Figures  10  and  11  is  referred  to  a line  perpendicular  to  the  crack 
rather  than  parallel  to  the  crack  as  in  Reference  16. 

The  anisotropy  factor  H increases  with  crack  length  for  the  CCT 
specimen  while  it  decreases  with  increasing  crack  length  for  the  DEC 
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TABLE  2.  ANISOTROPY  FACTOR  IN  CCT  ANGLE-PLY  SPECIMENS  OF 
T300/5208  GRAPHITE/EPOXY  [20] 


2a/W* 

or* 

0.1 

H 

0.2 

H 

0.3 

H 

0.4 

H 

0.5 

H 

0.6 

H 

B 

<°°>S 

1.000 

0.997 

0.991 

B 1 

H 

na 

(±10°) s 

1.001 

1.000 

0.997 

EES 

IE9 

<±15°>3 

1.002 

1.002 

1.003 

1.004 

1.005  ; 

(±20°) s 

1.003 

1.005 

1.010 

1.023 

1.030 

1.037 

(±25°)  g 

1.004 

1.010 

1.020 

1.033 

1.047 

1.061 

1.074 

(±30°) 8 

1.005 

1.016 

1.032 

1.052 

1.073 

.1 

1.113 

(±35°) s 

1.007 

1.021 

1.043 

1.070 

1.098 

1.145 

<±45“>. 

1.007 

1.025 

1.053 

1.087 

1.120 

1.152 

1.164 

<±55°>8 

1.002 

1.016 

1.037 

1.062 

1.087 

1.107 

1.149 

(±60°) s 

0.999 

1.008 

1.023 

1.040 

1.058 

1.080 

1.115 

(±65°) 8 

0.996 

1.002 

1.010 

1.020 

1.033 

1.073 

(±70°) b 

0.995 

0.997 

1.001 

1.006 

1.012 

1.021 

1.034 

(±75°) 8 

0.995 

0.995 

0.995 

0.995 

Esa 

0.997 

1.002 

(±80°) g 

EEa 

0.995 

0.992 

0.988 

n 

0.976 

(90°) s 

B 

0.995 

0.990 

0.983 

B 

0.964 

0.952 

*See  Figure  10 
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specimen.  This  suggests  that  the  effect  of  anisotropy  is  strongest 
when  the  crack  tip  is  near  a free  boundary.  For  both  specimens,  the 
severest  stress  intensity  occurs  when  the  plies  are  oriented  at  ±45°, 
which  agrees  with  the  behavior  observed  by  Synder  and  Cruse  [ 16]  in 
Figure  9.  The  lowest  stress  intensity  occurred  when  the  fibers  were 
either  aligned  with  the  crack  (a  “ 908)  or  perpendicular  to  the  crack  (a  = 0°)  . 
The  intensity  was  approximately  the  same  for  those  two  orientations. 
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2a/W* 

a* 


(°°)8 
(±l°°)a 
(±].50)a 
(±20°)g 
(±25°) g 
(±30°)a 
(±35°) g 
(±45°)g 
(±55°) g 
(±60°) 

<«•>: 

(±70°)g 

(±75°)a 

(±80°)8 

(90°) 


.945 
.969 
0.992 
1.019 
1.050 
1.083 
1.120 
1.142 
1.112 
1.092 
1.056 
1.022 
0.996 
0.973 
0.947 


*See  Figure  11. 


0.2 

H 

0.3 

H 

n 

ESS 

1.002 

1.026 

1.052 

1.055 

1.082 

1.086 

1.115 

1.112 

1.148 

1.135 

1.122 

1.116 

1.097 

1.089 

1.067 

1.059 

1.037 

1.029 

1.008 

1,005 

0.985 

0.986 

0.960 

0.967 

.972 
.987 
.003 
1.023 
1.049 
1.079 
1.107 
1.131 
1.098 
1.072 
1.045 
1.021 
1.002 
0.987 
0.973 


.979 

.991 

.003 

1.019 

1.041 


1.116 


1.035 

1.016 

1.001 

0.990 

0.981 


.986 

.994 

.002 

1.013 

1.029 

1.049 

1.069 

1.087 

1.064 

1.045 

1.027 

1.013 

1.002 


.991 
.996 
.001 
1.008 
1.018 
1.031 
1.043 
1.057 


1.012 

1.004 

0.998 

0.993 
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ANISOTROPY  FACTOR  (H) 


again  in  agreement  with  Reference  16,  Figure  8.  The  stress  intensity 
factor  for  a given  ply  orientation  angle  and  that  angle's  complement 
were  approximately  the  same.  For  example,  the  values  for  the  ±10°  and 
±80°  ply  laminates  were  almost  identical.  In  most  cases,  the  isotropic 
solution  underestimates  the  stress  intensity  factor.  Thus,  the  isotropic 
solution  tends  to  give  conservative  results  for  the  critical  stress  inten- 
sity factor  when  used  in  a fracture  toughness  test.  The  same  results, 
used  to  make  a residual  strength  calculation  for  a cracked  member  where 
the  toughness  is  known,  would  underestimate  the  stress  intensity  at  the 
crack  tip  and  thus  result  in  an  inflated,  nonconservative  value  for  the 
strength  of  the  component. 

The  dependency  of  the  stress  intensity  on  the  crack  tip-free  edge 
distance  was  illustrated  by  plotting  the  anisotropy  factor  H as  a func- 
tion of  the  nondimensionalized  free  edge  distance  £ for  both  the  CCT 
and  DEC  specimens.  The  free  edge  distance  is 

fcf  » W/2  - a 

for  the  CCT  specimen  and 


for  the  DEC  specimen.  The  results  for  some  representative  laminates 
are  shown  in  Figure  12. 

For  i^/W  > 0.3,  the  anisotropy  effect  is  virtually  the  same  for  the 

two  specimens.  H for  both  specimen  approaches  unity  at  approximately 
=*  0.45. 

The  stress  intensity  factor  depends  strongly  upon  the  shear  modulus 
The  axial  modulus  ratio  E^/E^  Play®  only  a secondary  role,  a 

fact  which  has  been  demonstrated  by  Synder  and  Cruse  [ 16]  as  well  as 
Konish  [20].  If  the  shear  modulus  for  an  angle  ply  laminate  is  accom- 
panied by  the  stress  intensity  factor,  it  increases  to  a maximum  at 
a * 45°.  Konish  demonstrated  this  for  various  2a/W  ratios  of  the  CCT 
and  DEC  specimens.  Because  results  are  similar  for  the  two  cases,  only 
the  CCT  specimen  is  shown  in  Figure  13.  The  relative  shear  modulus  G 
is  the  shear  modulus  for  the  angle  ply  laminate  divided  by_the  shear 
modulus  for  the  quasi-isotropic  (0/±45/90>s,  laminate.  A G greater  than 

unity  indicates  a laminate  relatively  stiff  in  shear,  while  a G less  than 
unity  indicates  a laminate  relatively  soft  in  shear.  It  can  be  seen  in 
Figure  13  that  for  laminates  stiff  in  shear  (G  > 1),  the  stress  inten- 
sity exceeds  the  isotropic  value;  for  laminates  soft  in  shear,  the 
stress  intensity  is  lower  than  the  isotropic  value. 
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LAMINATE 
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FREE-EDGE  DISTANCE  (Cf/W) 

Figure  12.  Anisotropy  factor  versus  free-edge  distance  for  represen- 
tative laminates.  Open  points  denote  center-cracked  specimens; 
filled  points  denote  DEC  specimens  [20]. 

Figure  14  shows  a comparison  of  results  obtained  by  Osias  and  Cruse 
[15]  and  Konish  [20]  for  the  CCT  specimen  with  fibers  at  ±45°.  The 
material  properties  were  slightly  different.  The  properties  of  Equations 
(31)  were  used  in  Reference  20  while  properties  of  Equations  (15)  were 
used  in  Reference  15.  The  Gl2  used  by  Osias  and  Cruse  [15]  was  approxi- 
mately twice  that  used  by  Konish  [20];  hence,  the  Osias-Cruse  results 
should  be  higher  than  those  of  Konish.  However,  that  is  not  the  case; 
the  Konish  results  are  slightly  higher.  There  is  little  practical 
difference  between  the  two:  compare  H " 1.10  to  H ■ 1.16  for  a/b  m 0.7. 
It  should  be  kept  in  mind  that  the  greatest  departure  from  the  isotropic 
value  occurs  for  the  present  case,  i.e.,  when  the  fibers  are  at  an  angle 
of  ±45°. 
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Figure  13.  Anisotropy  factor  versus  relative  shear  modulus  for 
center-cracked  specimens  [20], 

Thus  far,  the  discussion  has  indicated  that  the  anisotropy  effect 

on  the  stress  intensity  is  not  extremely  large.  It  must  be  pointed  out 

though  that  the  numerical  results  are  not  in  universal  agreement  with 

this.  For  example,  consider  the  boundary  collocation  results  of  Bowie 

and  Freese  [19]  for  the  square  orthotropic  CCT  specimen  in  Figure  15.  For 

values  of  E /E  greater  than  unity,  results  for  different  cracks  lengths 
x y 

are  only  moderately  lower  than  for  an  isotropic  case,  but  for  values  of 

E /E  less  than  unity,  that  is,  with  the  stiff  axis  perpendicular  to  the 
x y 

crack,  the  stress  intensity  increases  sharply  to  values  considerably  above 
the  isotropic.  Furthermore,  for  a/b  ■ 0.4,  the  results  have  been  veri- 
fied by  Atluri  et  al.  [18]  using  the  finite  element  technique.  The  reason 
for  the  behavior  exhibited  in  Figure  15  can  perhaps  be  explained  by  two 
effects.  First,  the  plate  is  square,  which,  compared  to  a long  plate, 
increases  the  finite  correction.  The  second  effect  is  related  to  the 
material  properties.  The  problem  was  formulated  in  terms  of  two  material 
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Figure  15.  Finite  correction  factor  for  square  plate  of 
orthotropic  material  [ 19] . 


r 
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constants,  and  p^»  which  are  combinations  of  the  four  orthotropic 

material  constants  E , E , G and  v , i.e., 

x’  y*  xy  xy’  * 


v*  -A71; 

h + h-J*  {y^r  + v2Giy  - v}172  • 


Elastic  constants  for  real  materials  were  not  chosen.  Instead,  relative 
values  of  P and  P„  were  used,  i.e.,  p “ 1 was  chosen,  which  leads  to 

P2- 

unspecified.  Their  values,  however,  are  implied  by  the  second  of 
Equations  (32).  In  fact,  if  a typical  value  for  Poisson's  ratio  is 
assumed  and  G is  calculated  in  terms  of  either  E or  E for 


E ^/E  . The  solution  was  then  completed  by  choosing  values  of 

rE  ) . This  left  the  shear  modulus,  G , and  Poisson' s ratio,  V , 
y7  xy  * xy* 


xy  x y 

E /E.  = 0.3,  G turns  out  to  be  larger  by  a factor  of  four  or  five 
x t xy 

than  a unidirectional  graphite  epoxy  laminate.  In  other  words,  picking 
specific  values  of  p^  and  p2  in  some  cases  implies  a shear  stiffness  which  is 


artifically  high  compared  to  real  materials.  It  has  already  been  noted,  based 
on  the  work  of  Konish  [20]  as  well  as  others  [ 15,  16]  , that  high  values  of 
shear  stiffness  lead  to  high  anisotropy  factors  as  shown  in  Figure  13. 


The  work  of  Konish  [20],  Snyder  and  Cruse  [16],  and  Osias  and  Cruse 
[ 15] , indicates  that  the  anisotry  factor  for  angle  plies  is  greatest  when 
a = ±45°.  A finite  element  solution  obtained  by  Mandell  [ 21]  disagrees  with 
this.  Figure  16  shows  a comparison  of  the  stress  intensity  from  the 
finite  element  technique  [21]  compared  to  the  boundary  integral  equation 
method  [20].  The  two  methods  agree  fairly  well  from  0°  to  approximately 
30°  and  from  75°  to  90°,  but  disagree  sharply  in  the  region  from  30°  to 
75°.  In  fact,  the  finite  element  results  show  the  stress  intensity  to 
decrease  to  a minimum  value  at  approximately  a ■ 50°,  near  the  place 
where  the  maximum  occurs  according  to  the  boundary  integral  equation 
method  [20] . The  reason  for  the  disagreement  is  not  known;  more  numeri- 
cal work  will  be  needed  to  settle  the  issue  completely.  It  is  believed 
that  because  the  boundary  integral  equation  technique  has  been  widely 
applied  to  anisotropic  crack  problems  by  several  researchers,  a great 
deal  of  confidence  can  be  placed  in  it. 

The  effect  of  plate,  length  to  width  aspect  ratio  for  a DEC  specimen 
was  illustrated  by  Atluri,  et  al.  [18].  The  results  for  fibers  either 
parallel  or  perpendicular  to  the  crack  are  shown  in  Figure  17  for  a 
length  to  width  ratio,  L/b,  of  4.  These  results  agree  with  the  behavior 
found  by  Konish  and  noted  already:  for  low  a/b  the  isotropic  solution 
would  overestimate  the  stress  intensity;  for  large  a/b  the  isotropic  and 
anisotropic  cases  converge.  Figure  18  for  L/b  ■ 1 (square  plate)  tells 
a different  story.  With  the  material's  stiff  axis  parallel  to  the  crack 
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FIBER  DIRECTION,  ±a  (deg) 

Figure  16.  K-Calibration  for  various  ply  orientations  (double-edge- 
notched  specimens,  graphite/epoxy,  2a/W  " 0.3). 


with  double  external  cracks.  Material  principal  directions  parallel 
to  geometric  axes  of  symmetry.  (L/b  =»  4)  [18]. 


the  anisotropic  stress  intensity  factor  is  as  much  as  507,  greater  than 
the  isotropic  value.  Whereas  the  isotropic  case  is  not  very  sensitive 
to  the  aspect  ratio,  a/b,  (compare  the  isotropic  solution  in  Figure  17 
with  that  in  Figure  18)  the  anisotropic  solution  is  very  sensitive.  While 
this  indicates  that  a great  deal  of  caution  must  be  exercised  when  using 
isotropic  solutions  for  unidirectional  laminates,  it  should  be  kept  in 
mind  that  L/b  = 1 is  a rather  severe  aspect  ratio— one  that  for  testing 
purposes  could  easily  be  avoided.  With  the  material's  stiff  axis  per- 
pendicular to  the  crack,  the  isotropic  solution  once  again  gives  a con- 
servative estimate  of  the  anisotropic  stress  intensity  factor. 

For  all  examples  considered  thus  far,  the  crack  line  has  coincided 
with  one  of  the  laminate's  principle  axes.  Atluri  et  al.  [18]  studied 
the  stress  intensity  behavior  in  the  DEC  specimen  with  the  principal 
axes  at  either  +45°  or  -45°  to  the  crack  line  as  shown  in  Figure  19. 
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Figure  18.  Stress  intensity  factor  in  an  orthotropic  tension  plate 
with  double  external  cracks.  Material  principal  directions  parallel 
to  geometric  axes  of  symmetry.  (iyb  = 1)  [18]. 

The  laminate  corresponds  to  a graphite-epoxy  material.  For  this  problem, 
which  lacks  tne  symmetry  of  the  previous  examples,  a value  of  exists 

also.  The  greatest  departure  of  from  the  isotropic  case  and  the 

largest  value  of  occur  at  approximately  the  same  location,  a/b  “ 0.4. 

It  should  be  noted  that  even  though  this  problem  represents  a severe 
example  in  terms  of  material  directionality  and  fiber  orientation  to  the 
crack,  for  a/b  > 0.6,  the  isotropic  and  anistropic  solutions  for  1C.  agree 
quite  closely. 

2.3  Mixed  Orthotropic  Laminates 

Results  so  far  have  been  restricted  to  either  unidirectional  or 
angle  ply  laminates.  In  practice,  a mixed  type  of  orthotropic  laminate 
is  frequently  encountered.  It  is  reasonable  to  expect  that  the  effect 
of  anisotropy  would  be  less  for  these  materials  because  they  are  less 
directional.  Figure  20  [15]  confirms  this,  where  a (±45)  laminate  for 
a CCT  specimen  is  compared  with  three  mixed  laminates.  The  isotropic 
solution  is  within  approximately  47,  of  the  three  mixed  laminates.  Data 
[ 20]  for  several  mixed  orthotropic  laminates  of  various  component  materials 
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Figure  19.  Stress  intensity  factor  in  an  orthotropic  tension  plate 
with  double  external  cracks.  Material  principal  directions  45°  to 
geometric  axes  of  symmetry  [ 18] . 

are  shown  in  Tables  4 and  5.  Again,  the  isotropic  solution  is  seen  to 
be  within  4%  of  correct  results  and  indeed  for  most  cases  within  17.. 

A number  of  problems  were  solved  by  Cruse  and  Osias  [ 15]  for  circular 
and  elliptical  holes  with  symmetric  radial  cracks  contained  in  plates  of 
mixed  laminates.  An  example  is  shown  in  Figure  21.  The  anisotropy  fac- 
tor for  a ±45°  laminate  is  compared  with  three  mixed  laminates.  Except 
for  small  crack  length,  a,  the  quantity  H is  near  unity  for  the  three 
mixed  laminates.  Again,  it  can  be  seen  that  the  (±45)  laminate  deviates 
the  most  from  isotropic  behavior,  being  approximately  107.  higher  than 
the  isotropic  solution.  It  is  important  to  note  the  behavior  for  the  short 
cracks.  As  the  crack  first  penetrates  the  hole's  edge,  the  anisotropic 
solution  either  overshoots  or  undershoots  the  isotropic  value  on  the  order  of 
107..  As  the  crack  grows  in  length,  the  solution  (for  the  mixed  laminates) 
approaches  the  isotropic  solution.  Thus,  Figure  21  indicates  that 
when  making  residual  strength  caluclations  for  short  edge  cracks,  a 
great  deal  of  caution  is  required.  If  most  of  the  fibers  are  perpen- 
dicular to  the  crack,  then  the  calculated  Isotropic  stress  Intensity 
should  perhaps  be  increased  by  approximately  107.  for  a conservative 
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TABLE  4.  ANISOTROPY  FACTOR  IN  DEC  SPECIMENS  OF  MIXED  ORTHOTROPIC 
LAMINATES  [20] 


a 

Material 

2a/W* 

(0°/±45o)g 

HTS 

H 

<0- /*«•>. 

Boron/Epoxy 

H 

(0°2/±45°)s 

HTS 

H 

(0°2/±45°)8 

Glass/Epoxy 

H 

• 

/"S  1 

e 

m 

• 

o 

s-/ 

0.1 

1.040 

1.010 

1.003 

MEM 

0.2 

1.037 

1.012 

1.012 

1.004 

MEM 

EaH 

1.014 

1.014 

1.006 

MEM 

BB 

1.013 

1.013 

1.006 

MM 

B 

1.014 

1.014 

1.008 

1.001 

1.018 

1.007 

1.006 

1.003 

WSSM 

1.011 

1.004 

1.004 

1.002 

warn 

See  Figure  11.  j 


f. 


TABLE  5.  ANISOTROPY  FACTOR  IN  CCT  SPECIMENS  OF  MIXED  ORTHOTROPIC 
LAMINATES  [20] 


a 

Material 

2a/W 

(0°/±45°)8 

HTS 

H 

(0°2/±^5°)8 

Boron/Epoxy 

H 

(0°2/±45°)8 

HTS 

H 

(0°2/±45°)s 

Glass/Epoxy 

H 

(0°4/±45°) 

HTS 

H 

0.1 

1.002 

1.002 

1.001 

mm 

mm 

1.006 

I 

1.003 

1.001 

19 

1.012 

1.005 

1.005 

■E9 

0.4 

1.008 

1.008 

■S3 

0.5 

1.011 

1.011 

1.004 

0.6 

1.037 

1.014 

1.014 

1.005 

■S3 

0.7 

1.045 

1.017 

1.017 

1.006 

■9 

*See  Figure  10. 
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Figure  21.  Anisotropic  effects  for  small  circular  notched  plate  t 15] 


strength  estimate.  For  engineering  calculations  a similar  approximate 
correction  could  also  be  applied  to  the  isotropic  solution  for  (£45) 
laminates  (except  for  extremely  short  edge  cracks) , but  for  most  mixed 
laminates  the  isotropic  solution  appears  to  be  sufficiently  accurate 
for  most  practical  applications. 

The  effect  of  short  cracks  at  the  edge  of  a hole  is  also  illustrated 
in  Figure  22  for  a specific  hole  radius  of  0.15  in.  in  several  mixed 
laminates.  The  values  were  calculated  by  Bowie  and  Freese  but  were 
included  via  personal  communication  in  a report  by  Olster  and  Woodbury 
[ 22] . The  results  [ 22]  were  used  for  making  residual  strength  calcula- 
tions of  panels  damaged  by  ballistic  impact  (Chapter  5) . For  all  lami- 
nates, the  material  properties  increase  the  stress  Intensity  above  the 
isotropic  value;  but  as  the  crack  grows,  the  laminate  values  approach 
the  isotropic  solution. 
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Chapter  3.  FUNDAMENTAL  FRACTURE  CRITERIA 
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As  seen  in  Chapter  2,  within  the  framework  of  fiber-reinforced 
laminate  theory  where  it  is  assumed  that  filamentary  laminates  behave 
in  a homogeneous  and  anisotropic  manner,  the  stress  analysis  for  cracks 
presents  no  serious  obstacle.  The  physical  process  of  crack  extension 
in  a laminate,  however,  tends  to  be  very  complicated;  many  modes  of  fail- 
ure are  possible.  In  a unidirectional  composite,  the  crack  most  frequently 
extends  collinearly  in  the  direction  of  the  fibers;  in  a cross-ply  laminate, 
the  crack  frequently  extends  in  an  irregular  manner  following  the  fibers 
of  one  ply  for  a distance  and  then  turning  to  follow  the  fibers  of  a cross 
ply.  Delamination  of  the  lamina  may  occur  in  the  crack-tip  region. 

Resulting  fracture  surfaces  range  from  smooth  to  the  "shaving  brush"  type. 
Obviously,  the  heterogeneity  of  the  laminate  strongly  influences  the  crack 
behavior.  In  light  of  this,  the  question  is:  Is  it  possible  to  character- 
ize the  fracture  behavior  of  filamentary  composites  using  the  theory  of 
linear  elastic  fracture  mechanics  (LEFM)  along  with  its  assumption  of 
homogeneity?  Considerable  supporting  evidence  of  this  idea  now  exists  — 
evidence  shows  that  conventional  fracture  mechanics  can  indeed  be 
used  to  predict  the  strength  of  cracked  (notched)  laminates. 

3.1  Fracture  of  Unidirectional  Composites 

The  discussion  in  this  section  is  limited  to  the  fracture  behavior 
of  a unidirectional  laminate  containing  a crack  oriented  along  the  fiber 
direction.  Admittedly  a special  case,  the  validity  of  LEFM  applied  to 
unidirectional  laminates  is  fundamental  to  the  more  general  question  of 
the  validity  of  LEFM  for  built-up  laminates.  A series  of  tests  was  con- 
ducted by  Uu  [5,  6]  to  determine  under  what  conditions  LEFM  could  be 
applied  to  the  fracture  of  unidirectional  composites.  His  specific 
experimental  objectives  were  as  follows; 

a)  To  determine  the  mode  of  crack  extension  in  the  presense  of 
both  symmetric  and  skew  symmetric  loads. 

b)  To  determine  if  the  concept  of  critical  kj  and  k_  could  be 
applied. 

c)  To  determine  a fracture  criterion  for  the  case  of  combined  mode. 

Two  materials  were  used  in  the  tests:  balsa  wood  sheets  and  fiber  glass 
reinforced  epoxy  sheets  known  as  Scotch  Ply  1002. 

The  balsa  wood  sheets,  0.0625  in.  thick,  and  the  Scotchply  sheets, 
0.05  in.  thick,  were  both  prepared  with  artificial  flaws  oriented  along 
the  grain  or  fiber  direction.  The  balsa  wood  plates  were  loaded  to  fail- 
ure using  the  four  load  paths  shown  in  Figure  23;  and  the  Scotchply  was 
loaded  to  failure  using  load  paths  Mo.  1,  2,  and  4.  Thus,  the  flawed 
plates  were  tested  with  tension  normal  to  the  crack;  a combination  of 


Figure  23.  Load  paths  used  by  Wu  [5,  6]. 

tension  and  shear,  and  pure  shear.  To  prevent  buckling  in  shear,  the 
balsa  wood  sheets  were  clamped  between  glass  plates;  the  Scotchply 
sheets  were  clamped  between  Teflon-lined  plastic  plates.  There  were 
several  reasons  for  oreinting  the  flaw  along  the  grain.  The 
orientation  of  the  flaw  must  be  known  to  find  k^  or  k^.  and  the 

preferred  orientation  of  a natural  flaw  is  along  the  fibers.  Also  there 
was  a consideration  of  displacements;  the  stress  intensity  factor  for 
the  anisotropic  case  is  consistent  with  the  isotropic  case  insofar  as 
stress  distribution  is  concerned,  but  not  for  displacements,  i.e.,  a 
symmetric  load  can  produce  skew-symmetric  displacements  and  a skew  sym- 
metric load  can  cause  crack  opening.  This  effect  disappears  if  the 
crack  is  along  one  of  the  material's  principal  directions.  Moreover, 
the  extreme  of  the  critical  stress  intensity  factor  should  occur  when 
the  crack  is  oriented  along  one  of  the  principal  directions. 

For  all  load  paths  Including  pure  shear,  the  crack  propagated 
collinearly  along  the  direction  of  the  fibers.  To  determine  if  critical 
values  of  k^  and  k^  are  indeed  material  constants,  it  is  necessary  to 

determine  their  form  in  terms  of  crack  length  and  load.  For  example,  for 

load  path  No.  1,  ■ 0y/~&,  ■ 0 (recall  that  m Jii  Kj  and  ■ 

yir  k^).  Thus,  if  critical  k^  is  a material  constant,  then  the  critical 

stress  should  be  proportional  to  (critical  crack  length)*^.  Thus,  a plot 
of  critical  stress  versus  critical  crack  length  on  a log-log  scale  should 
result  in  a straight  line  with  a slope  m of  -1/2. 
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Wu  plotted  the  critical  stress  (either  the  critical  normal  stress 
or  the  critical  shear  stress  component  as  appropriate  to  the  load  path) 
versus  critical  crack  length  and  fitted  the  data  with  a straight  line  of 
slope  m.  Figure  24  shows  the  results  fcr  load  path  No.  1 (tension)  for 
both  balsa  wood  and  Scotch  Ply  1002.  The  values  of  m for  all  load  paths 
ranged  from  -0.47  to  -0.55  for  balsa  wood  and  from  -0.46  to  -0.50  for 
Scotch  Ply,  in  both  cases  very  near  the  value  of  -1/2.  Thus,  it  was  con- 


1 .2  .3  .4  .5  .6  .7.8.9  1 .15  .2  .3  .4  .5  .6  .7  .8.9  1 

CRITICAL  HALF  CRACK  LENGTH,  a (in.)  CRITICAL  HALF  CRACK  LENGTH,  a (in.) 


Figure  24.  Result  of  tension  tests  for  balsa  wood  and  fiber  reinforced 
epoxy  (Scotchply  1002)  [5]. 

An  extending  crack  in  an  isotropic  plate,  typically  turns  to  a direc- 
tion normal  to  the  greatest  applied  principal  stress.  But  the  crack  in 
these  experiments  extended  collinearly  under  a combination  of  Mode  I and 
Mode  II.  Thus,  unlike  the  isotropic  case,  it  becomes  important  to  deter- 
mine a fracture  criterion  for  combined  mode  extension.  Wu  [6]  considered 
the  maximum  energy  release  rate  criterion  of  Irwin  [23]. 

*1  + *11  " 8C  " const*nt  <33> 
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where  from  Equations  (14) 


8II  " * 


. 2 /a22all  /*22  . 2 a12  + a66T 
.2  f Fri  . 2 a12  + a66l 

^[A:  2 *n  j 


Substituting  for  and  into  Equation  (33)  yields 


.2  . All  .2 
kT  + / ^ = c 
1 J a22  “ 


where  the  right-hand  side  has  been  redefined  as  C',  a constant.  If  Equa- 
tion (34)  is  true,  then  it  should  hold  for  both  of  the  following  cases: 


pure  tension,  = k^c  k^  = 0 

pure  shear,  k^  « 0 = kjIC 


These  two  conditions  in  Equation  (34)  result  in 


all  _ / kic  \ 
a22  YTIc) 


which,  when  substituted  into  Equation  (34),  yields 


fe)  •&)- 


This  is  the  energy  release  rate  criterion  in  terms  of  stress  intensity 
factors.  Figure  25  shows  the  values  of  failure  of  and  k^  for  the 

Scotchply  material  tested  by  Wu.  It  is  instructive  to  compare  the  data 
with  the  strain  energy  release  rate  criterion.  Equation  (33)  or  equiva- 


lently, Equation  (34).  Let  the  right-hand  side  be  given  by  k^  and  use 

1/2  1C 

the  value  reported  by  Wu  for  k^  of  930  psi-(in.)  . Further , Wu  reported 


an/*22  " 0.33.  These  values,  inserted  into  Equation  (34),  result  in 
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the  dashed  line  on  Figure  25. 


It  is  apparent  that  g^.  + ■ constant 


does  not  adequately  describe  the  fracture  behavior  in  the  presence  of 
combined  mode. 

For  purposes  of  fitting  the  data,  the  results  were  plotted 
in  nondimens ional  form  as  shown  in  Figure  26 , The  equation 


(36) 


fits  the  data  very  well  and  is  shown  as  the  solid  dark  line.  Equation  (25) 
is  represented  by  the  dotted  line.  While  Equation  (36)  provided  a good  fit  to 
the  data  for  both  balsa  wood  and  Scotchply,  it  should  be  kept  in  mind 
that  Equation  (36)  is  basically  only  an  empirical  equation  which  explains 
nothing  fundamental  about  the  fracture  mechanism.  On  the  other  hand,  the 
g approach  represented  by  Equation  (35)  is  fundamental.  Its  failure 
to  represent  the  data  must  be  considered. 

During  the  fracture  tests,  the  crack  extension  process  observed.  It 
was  noticed  that,  though  the  crack  basically  followed  the  fiber  direction, 
it  did  not  extend  in  a completely  planar  fashion,  that  it  intermittently 
skipped  back  and  forth  across  bundles  of  fibers  as  shown  in  Figure  27. 

This  resulted  in  the  upper  and  lower  fracture  surfaces  being  connected 
by  fibers.  Corten  [12]  discussed  the  idea  that  these  fibers  acted  as 
shear  connectors  for  Mode  II  displacements  but  that  because  the  fibers 
were  parallel  to  the  fracture  surface,  they  carried  almost  no  load  nor- 
mal to  the  crack  surface.  The  effect  of  the  shear  connectors  was  to  cause 
an  apparent  value  of  k^^,  larger  than  the  actual  k^^,. 

One  other  point  about  the  g approach.  Equation  (35),  must  be  dis- 
cussed. Lauratis  [24]  conducted  some  off-axis  tension  tests  on  Scotchply 
which  contained  only  natural  flaws  in  the  form  of  bubbles,  voids 
fiber-matrix  delamination  from  lack  of  wetting,  etc.  She  obtained  good 
correlation  of  Equation  (35)  with  her  test  results.  The  question  is,  if 
Equation  (35)  would  not  correlate  forWu's  experiments,  did  it  correlate  for 
Lauraitis?  Unfortunately,  Lauraitis  failed  to  discuss  this  point  even  though 
she  used  Wu's  data  for  k^,  and  a^/a 22*  ®ne  can  8Peculate  38  to  why  she  found 

good  correlation.  The  answer  may  lie  in  the  model  she  chose.  She  modeled 
the  microflaws  as  tunnel  cracks  and  inferred  the  microflaw  crack  length 
from  her  tension  tests.  A penny-shaped  flaw,  (as  she  mentions)  or 
an  elliptical  flaw  may  well  have  been  more  representative  of  the  actual  flaw 
geometry.  A different  flaw  model  would  have  changed  not  only  the  magnitude 
of  but  would  also  have  changed  the  relative  amounts  of  k^  and  k^. 

Thus,  the  tunnel  shaped  model  may  very  well  have  been  fortuitious  and 
a more  realistic  crack  shape  might  have  actually  given  poorer  corre- 
lation. In  any  case,  the  answer  is  not  clear  at  this  time  and  more 
work  will  be  needed  to  settle  the  question.  In  the  course  of  her  work 
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Lauraitis  [24]  established  the  existence  of  mlcroflaws  oriented  along 
the  fiber  direction  in  the  fiber  glass  epoxy. 

The  general  formulation  of  the  combined  mode  problem  will  be  dis- 
cussed further  later  on;  the  point  to  be  made  here  is  that  Wu  estab- 
lished the  validity  of  applying  fracture  mechanics  to  orthotropic  plates, 
more  specifically,  to  unidirectional  fiber  composites  with  fiber-aligned 
cracks.  Others  have  confirmed  this.  Konish  et  al.  [ 7,  8]  tested  graphite- 
epoxy  using  three  point  bending,  Slepetz  and  Carlson  [25]  tested 
glass  epoxy  and  two  graphite  epoxies  using  the  compact  tension  specimen. 

For  the  case  where  the  initial  crack  is  aligned  with  unidirectional 
fibers,  the  experimental  data  from  tests  on  several  fiber  composite 
systems  agree:  the  fracture  path  extends  generally  along  the  fiber 
direction,  the  critical  stress  intensity  factors  are  material  constants, 
and  in  spite  of  the  heterogeneity  of  the  material,  LEFM  can  be  usefully 
applied  to  predict  the  fracture  of  such  systems. 

3.2  Fracture  of  Multidirectional  Laminates 

This  section  deals  with  the  fracture  of  multidirectional  laminates 
of  varying  stacking  sequences  and  angles.  Some  information  on  fracture 
of  unidirectional  composites  where  the  original  crack  plane  does  not 
coincide  with  the  fiber  direction  will  also  be  included.  Because  of  the 
many  possible  lay-up  patterns  for  laminates,  fracture  surfaces  are  varied 
and  frequently  complex  in  shape.  General  statements  are  more  difficult 
to  make  than  in  the  last  section;  results  for  several  laminates  of 
different  stacking  sequences  and  angles  must  be  considered.  In  doing 
this,  there  will  be  less  emphasis  on  fracture  toughness  values  (Chap- 
ter 4)  than  on  the  degree  of  applicability  of  LEFM  to  such  systems. 

• 

The  applicability  of  LEFM  to  composites  has  been  investigated  by 
Konish  and  Cruse  [8]  and  Konish  et  al.  [7,  9].  They  ran  a number  of 
experiments  on  cracked  specimens  of  different  configurations.  The 
specific  purposes  were  to  determine  whether  the  fracture  path  was  deter- 
mined by  the  geometry  of  the  initial  crack  and  the  loading  or  by  the 
material  orientation  and  to  investigate  whether  LEFM  could  be  usefully 
applied  to  composites. 

In  keeping  with  the  objective,  the  test  procedures  followed  those 
developed  within  the  framework  of  conventional  fracture  mechanics.  They 
employed  the  three-point  bend  specimen  shown  in  Figure  28.  The 
material  was  a graphite  epoxy.  To  establish  reproducibility  five  tests 
each  were  run  on  a unidirectional  laminate  (a  ■ 0°)  and  a multidirectional 
laminate  [a  = (,0°/±A5°/90°)  b]  . The  crack  length,  a,  for  both  materials 

was  0.4  in.  Single  tests  were  then  run  for  a m 0°,  45°,  90°  (±45) 

and  (0o/±45o/90o)8, all  with  starter  cracks  of  0.2,  0.4,  and  0.6  in.  The 

short  crack  lengths  were  included  to  permit  evidence  of  material  dominance 
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Figure  28.  Three-point  bend  specimen  showing  fiber  direction. 


to  develop.  The  beams  were  loaded  to  failure  and  the  critical  stress 


intensity  factors,  called  here,  were  computed  for  all  three  crack 


lengths.  The  values  were  determined  from  the  isotropic  formulas. 


Results  for  along  with  g^,  the  critical  strain  energy  release  rate. 


are  shown  in  Table  6.  To  be  a material  constant,  should  be  independ- 


ent of  the  crack  length.  While  perhaps  not  conclusive,  results  indicate 


that  within  experimental  error  is  essentially  constant  for  each  lami- 


nate. It  should  be  noted  that  for  the  45°  laminate,  the  presence  of 
tension-shear  coupling  would  induce  both  Mode  1 and  Mode  II  behavior  at 


the  crack-tip.  will  reflect  the  presence  of  both  modes  although  the 


relative  amount  of  Kj.  and  was  not  known.  The  grouping  of  the  tough- 

ness values;  the  toughness  of  the  80°  and  45°  laminates,  is  an  order 
of  magnitude  less  than  the  other  three  because  the  90°  and  45®  lamin- 
ates fractured  through  the  matrix  whereas  the  others  fractured  across 
fibers. 


Figure  29  shows  the  load-displacement  curves  for  the  reproducibility 
tests  for  the  a - (0°/±45°/90°)  laminate.  As  a matter  of  convenience, 


the  cross-head  displacement  was  monltered  rather  than  the  crack-opening 
displacement  as  specified  by  ASTM  for  the  testing  of  metals  [26].  The 
shape  of  the  curves  in  Figure  29  is  important.  The  curves  exhibit  linear 
behavior  over  most  of  their  range  with  only  a very  slight  amount  of  non- 
linearity preceding  the  maximum  load  or  instability.  In  the  testing  of 
composites,  there  is  a crack-tip  damage  zone  which  is  somewhat  analogous 
to  the  crack-tip  plastic  zone  for  metals.  This  dasiage  zone  is  due  to 
such  things  as  splitting  parallel  to  the  fibers  and  delamination  between 
plies.  A great  amount  of  such  crack-tip  damage  causes  increasing  com- 
pliance with  increasing  load,  leading  to  a highly  nonlinear  load- 
displacement  curve.  In  the  extreme  case,  the  fracture  would  be  a slow 
tearing  type  and  the  application  of  brittle  fracture  mechanics  becomes 
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TABLE  6.  CRITICAL  STRESS  INTENSITY  FACTORS  FOR  SEVERAL 
GRAPHITE  EPCKIES  [7] 


Fiber 

Orientation 
Angle,  Q 

V 

lb/in.2 Tin.  X 10"3 

Kq,  lb/in.2 

yin.  X 10-3 
(average) 

8q» 

lb/in.2 

a ■ 0.2  in 

. a ■ 0.4  in.  a m 0.6  in. 

0° 

1 

28.8 

36.3 

32  +m 

32  ’ 6 -117. 

■ 

90° 

1.66 

1.46 

2 

+6.37. 
1,56  -6.37. 

■ 

45° 

0.6903 

2.22 

2.39 

2 30  +3*^/o 
-3.87. 

HH 

(±45°) 8 

18.5 

18.5 

16.3 

+4.87. 
1/>7  -9.47. 

45.0 

(0°/±45o/90°)a 

23.5 

21.7 

20.5 

,,  „ +7.37. 
-8.67. 

55.1 

1 

Notes: 

1.  Specimen  was  crushed  before  crack  propagation  occurred. 

2.  Instrumentation  failure. 

3.  This  value  omitted  when  calculating  ave. 

4.  No  gq  available  because  the  crack  propagated  in  a mixed  mode,  which 
could  not  be  directly  uncoupled. 


suspect.  The  curves  in  Figure  29,  however,  exhibit  only  a slight  amount 
of  nonlinearity,  indicating  that  the  failures  are  essentially  brittle, 
i.e.,  catastropic  with  little  prior  warning,  and  that  the  usual  ideas 
of  brittle  fracture  mechanics  apply. 

The  tests  [7]  were  run  for  several  laminate  types  and  the  fracture  sur- 
faces varied.  In  the  absense  of  photographs,  a brief  description  is 
given  of  each  type  as  follows: 

a)  a ■ 0° — The  path  of  crack  growth,  though  somewhat  irregular,  was 
roughly  coplanar. 

b)  Of ■ 45° —The  crack  grew  along  a plane  containing  no  fibers, 
i.e.,  at  45°,  to  the  starter  crack,  and  the  resulting  fracture  surface 
was  very  nearly  planar.  Crack  propagation  took  place  here  under  a com- 
bination of  Kj  and  Kjj. 


. 
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Figure  29.  Load  versus  cross-head  displacement  for  five  three-point 
bend  specimens  (reproducibility  tests)  of  a ( 0“/±45®/90° ) graphite 
epoxy  laminate  (a  ■ 0.4  in.)  [7].  * 


c)  a ■ 90*  —The  crack  grew  along  a plane  containing  no  fibers 
coplanar  with  the  starter  crack. 

d)  a ■ (±45®)^  — For  a ■ 0.2  in.,  the  crack  almost  immediately 
turned  45*  from  the  started  direction.  For  a ■ 0.4  in.,  there  was  a 


greater  tendency  for  coplanar  growth.  For  a - 0.6  in.,  crack  dominance 
was  apparent  with  almost  coplanar  growth.  For  both  a ■ 0.2  and  0.4, 
zig-zagging  was  apparent. 

e)  a ■ (0°/±45o/90°)g —Crack  growth  was  coplanar. 

That  crack  dominance  of  the  fracture  process  decreases  for  the 
shorter  crack  lengths  is  indeed  apparent  from  the  preceding  fracture 
surface  descriptions.  This  leads  to  questions  about  minimum  crack  size 
for  which  the  criterion  will  apply.  Descriptions  of  near  microscopic 

details  of  crack  growth  in  graphite  epoxy  and  glass  epoxy  laminates  have 

been  given  by  Mandell  [21],  In  that  study,  a minimum  crack  size  of  0.1  in. 

for  two  laminates:  (90/0/90/5)  and  (45/0/-45/90)  was  given  above  which 

8 8 

fracture  was  expected  to  be  controlled  by  the  criterion  and  below  which 

the  application  of  fracture  mechanics  becomes  doubtful.  The  minimum 
initial  crack  size  for  which  the  crack-tip  stress  field  will  control  the 
fracture  depends  upon  the  lay-up  of  the  particular  laminate  and  strickly 
speaking  would  require  determination  for  each  case.  As  an  alternative, 
the  minimum  crack  size  [21]  might  be  adopted  for  similar  laminates  as  an 
approximate  value  useful  in  design  situations. 


3.3  A Lamination  Theory  of  Fracture  [15,  27] 

By  properly  controlling  the  orientation  of  the  constituent  plies, 
it  is  possible  to  build  a laminate  in  such  a way  as  to  produce  predict- 
able stiffness  properties,  that  is,  to  tailor  a laminate  from  its  constituent 
plies  for  a specific  application.  It  would  thus  be  advantageous  to  be 
able  to  predict  the  fracture  toughness  or  the  critical  strain  energy 

release  rate  g^  for  the  laminate  from  the  values  of  the  critical  energy 

release  rates  for  the  constituent  angle  plies,  g^.  A relationship  for 

doing  this,  developed  by  Cruse  [15,  27],  shows  that  the  energy  released 
by  coplanar  crack  growth  at  fracture  is  the  sum  of  the  energies  released 
by  each  angle  ply  constituent  at  fracture.  The  relationship  is  restricted 
to  laminates  of  angle  ply  constituents  and  is  based  on  the  assumption  of 
coplanar  crack  growth,  which  requires  fracture  of  the  fibers.'  As  noted 
in  the  previous  article,  the  energy  required  to  fracture  a simple  ply 
by  matrix  failure  is  much  less  than  for  an  angle  ply  material. 

Tt\e  crack  is  modeled  as  a narrow  ellipse  (Figure  30)  with  a/b  « 1. 
The  stress  at  the  ends  of  the  ellipse  (x^  ■ ±a)  is  given  by 


_2  . a 
a 1 b 


[VlJ 

L 41M2  J 


where  |i^  and  ^ are  the  roots  of  the  characteristic  Equation  (17).  This 
equation  is  rewritten  as 
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Figure  30.  Elliptical  crack  model  [27], 
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and  the  last  term  is  neglected  in  comparison  with  Equation  (24)  is 
now  recalled  for  g , the  strain  energy  rate.  Substituting  Equation  (38) 
into  Equation  (24)  results  in 


8 - T a22  \k)  - w 

where  the  subscript  I has  been  dropped  from  g.  Assuming  that  the  stress 
intensity  for  an  anisotropic  specimen  is  the  same  as  for  an  isotropic 
specimen  of  identical  geometry,  the  following  is  written: 

Kj  - a Y Jl  (40) 

where  Y is  the  finite  correction  factor.  Substituting  Equation  (40)  into 
Equation  (39),  the  following  is  obtained: 

g - ^ Y2  a b (a22  a2)  . (41) 

At  the  tip  of  the  ellipse,  the  state  of  stress  is  uniaxial,  i.e.,  ■ 

t12  ■ 0.  Since  a22  ■ l/E22,  then  a22  a2  iS  th*  local  hoop  atrain  e2‘ 


In  the  limit  for  vanishing  b,  e2  is  unbounded  but  the  product  be2  is 
bounded  [27].  Defining  be2  to  be  5;  Equation  (41)  becomes 

g - Y2  a 6/2  . (/ 


The  quanity  6 is  interpreted  by  Cruse  as  a crack-tip  stretching,  physically 
related  to  the  fracture  process.  A fractographic  examination  by  Cruse 
and  Stout  [28]  indicated  that  for  a graphite  epoxy  laminate,  the  values 
of  &R  for  the  constituent  angle  plies  were  essentially  the  same.  This 


means  that  there  is  no  significant  crack  growth  in  any  ply  prior  to  frac- 
ture. Thus,  one  can  set  8^  ■ 5.  Furthermore,  let  the  value  of  6 at 

fracture  be  8^  and  the  stress  at  fracture  be  Op.  The  stress  trR  on  the 

constituent  angle  ply  of  thickness  is  related  to  the  overall  stress 

a on  the  laminate  thickness  B by 


N 

*B  " I <*KhK  (43) 

k«l  K K 


where  N is  the  number  of  constituent  angle  plies.  Substituting  for  a 
from  Equation  (43)  into  Equation  (42)  and  letting  fracture  conditions  apply 
results  in 


(44) 


Similai^y,  the  fracture  criterion  for  a constituent  angle  ply  can  be 
written: 


- A*  V2 


(45) 


Using  the  preceding  observation  that  8 

CJK. 

bined  with  Equation  (44)  to  yield 


Equation  (45)  can  be  com- 


(46) 


This  shows  that  the  energy  released  by  coplanar  crack  growth  is  the 
sum  of  the  energies  released  by  the  angle  ply  constituents,  a result 
which  perhaps  is  intuitively  apparent.  It  must,  however,  be  kept  in 
mind  that  Equation  (46)  applies  only  to  Mode  I,  i.e.,  the  crack  propa- 
gates in  its  plane,  and  the  damage  area  at  the  crack  tip  due  to  splitting 
between  fibers  and  delamination  between  plies  is  small.  Cruse  [27] 
applied  his  results,  Equation  (46)  to  the  data  from  Reference  7,  which 
are  shown  in  Table  6.  The  value  of  gp  for  the  (0/±45/90)  laminate  is 

computed  from  Equation  (46)  to  be: 

gp  - 117  (|)  + 45  (i)  + 0.9  (£)  + 52  in. -lb/in. 2 
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Compared  with  55.1  in. -lb/in.  , the  value  is  low  by  only  67..  Consider** 
ing  that  in  a laminate  failure,  there  are  probably  other  energy  consum- 
ing damage  processes  not  included  in  Equation  (46)  , the  results  are 
remarkably  close.  Cruse  and  Osias  [ 15]  have  checked  the  model  against 
the  gq  results  from  three-point  bending  for  a number  of  other  laminates, 

including  (90/±45),  , (0/±45) , and  (O^^S)  . The  g^  predicted 


from  Equation  (46)  was  for  most  cases  within  57.  of  the  measured  g^.  Equal 

accuracy  for  many  other  laminates  would  not  be  expected.  As  Cruse  [27] 
points  out,  if  the  crack  grows  by  delamination  followed  by  independent 
ply  failure,  then  Equation  (46)  would  not  apply. 


3.4  The  Energy  Density  Theory  [29,  30] 

As  already  noted,  the  Grif fith-Irwin  equation  which  relates  and 

kj  is  based  upon  collinear  crack  propagation.  Further,  the  idea  of  total 

strain  energy  release,  g^  + g^  = gc>  for  combined  mode  problems  depends 

upon  collinear  crack  propagation,  which  as  already  shown,  is  often 
contrary  to  experimental  evidence.  In  an  effort  to  circumvent  this  diffi- 
culty, Sih  and  Chen  [30,  31,  32]  have  devised  a fracture  criterion  that 
constitutes  a fundamental  departure  from  the  classical  theory  of  fracture 
mechanics.  The  method , which  employs  the  familiar  quantities  kj.  and  k^, 

is  based  upon  the  strain  energy  density  outside  a small  core  region 
surrounding  the  crack  tip.  With  this  method,  it  is  possible  to  calculate 
the  initial  crack  propagation  direction  as  well  as  the  critical  load. 

For  composites,  the  method  applies  to  unidirectional  laminates;  it  assumes 
that  fracture  occurs  by  crack  propagation  in  the  matrix. 

The  energy  density  theory  must  first  be  explained  before  its  appli- 
cation to  composites.  A basic  idea  of  the  theory  is  that  the  continuum 
mechanics  solution,  on  which  k^.  and  k^  are  based,  stops  short  of  the 

crack  tip  by  a distance  r , Figure  31. 


rQ  - RADIUS  OF  CORE  REGION 

Figure  31.  Core  region  surrounding  the  crack  tip. 
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The  radius  rQ  is  approximately  the  same  order  of  magnitude  as  the  crack- 

tip  radius.  Let  the  crack  tip  be  embedded  in  an  isotropic  (matrix)  mate- 
rial having  Young' 8 modulus,  E,  Poisson's  ratio,  v,  and  shear  modulus, 

G ■ E/[2(l  + V)] . Now  the  strain  energy  dw/dA  stored  by  the  element 
A A ■ rA0Ar  is  considered: 


dw  _ 1 / 

dA  2E  V 
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where  the  stress  components  r and  t are  zero  for  plane  strain.  The 

xz  yz 

isotropic  crack-tip  stresses  for  Mode  I loading  are  [1,  3], 
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polar  coordinates.  Upon  substituting  the  preceding  stresses  into  Equa- 
tion (47) , the  energy  density  can  be  written  as 


62 


dA  “ + 2 ai2kXklI  + 


The  coefficients  a^,  and  for  plane  strain  are  given  by 


ail  " 16G^3  ” 4v  “ coa  + co® 


°12  " lfe(2  8ln  6)  lcos  9 “ U " 2v)] 


a22  " T6G^4^  “ vX1  “ cos  e>  + (1  + cos  0)  (3  cos  0-1)] 


The  local  strain  energy  density  function  dw/dA  is  proportional  to  ; 

as  r becomes  smaller  and  smaller  dW/dA  becomes  larger  and  larger.  The 
intensity  of  the  strain  energy  is  referred  to  as  the  strain  energy  den- 
sity factor,  denoted  by  S,  where 

S " “ll1^  “ 2 0£l2kIkII  “ a22kII  * (52) 

This  is  the  key  quantity  in  the  application  of  the  theory.  It  depends 
upon  0;  thus,  it  is  direction  sensitive  and  its  role  in  characterizing 
the  strain  energy  is  somewhat  analogous  to  the  role  of  or  k^  in 

characterizing  the  stress  field.  In  the  energy-density  theory,  the  two 
hypotheses  of  unstable  crack  propagation  are  as  follows: 

a)  Crack  initiation  takes  place  in  the  direction  of  the  stationary 
value  of  the  strain-energy  density  factor,  that  is, 


§ - °.  at  0 - 0c 


and  the  crack  extends  in  the  direction  0 . 

o 

b)  Crack  extension  occurs  when  the  strain-energy  density  factor 
reaches  a critical  value,  i.e.,  when 

S(k1,kII)-Sc  , for  0 - 0o  . (54) 

The  difference  between  S and  S£  must  be  clearly  understood.  For  a given 
value  of  0,  S can  be  calculated  in  terms  of  crack  length  and  applied  load 
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through  kj  and  k^j.  On  the  other  hand,  is  a characteristic  of  the 

material,  a measure  of  the  material's  resistance  to  fracture.  The  advan- 
tages of  the  S-theory  are  that  it  is  applicable  in  the  presence  of  both 
Mode  I and  Mode  II  loading  and  it  will  predict  the  initial  direction  of 
crack  propagation. 


Sih  and  Chen  [30]  applied  the  theory  to  Wu's  data  [5,  6]  for  glass 
epoxy.  They  adopted  the  model  shown  in  Figure  32.  A crack  length  of 
2a  is  embedded  in  the  matrix  material  parallel  to  the  fibers,  where  the 
clear  distance  between  the  fibers  is  2h.  For  finding  the  stress  intensity 


factors,  they  assumed  that  the  material  at  a distance  greater  than  h from 
the  crack  is  homogeneous  and  orthotropic,  with  properties  E , E,.,  v19 

1 *m  lfc  f 


V21  and  G12*  Parallel  and  perpendicular  to  the  crack.  Figure  32.  The 


crack  is  contained  in  an  isotropic  layer  of  thickness  2h  between  these 
two  orthotropic  materials  as  shown  in  Figure  33. 


a 


Figure  32.  Fiber  composite  in  tension  with  crack  parallel  to  the  fibers. 

The  average  clear  spacing  2h  between  the  fibers  can  be  calculated  in 
terms  of  the  fiber  volume  fraction  v^  by  assuming  that  the  fibers  are 

distributed  uniformly  throughout  the  matrix.  Consider  the  fiber  of 
radius  R shown  in  Figure  34  surrounded  by  its. portion  of  the  matrix 
material.  The  dotted  lines  represent  the  centerlines  between  fibers. 

For  a unit  thickness  perpendicular  to  the  paper,  the  total  volume  VT  is 

VT  - (h  + 2R)2  (55) 
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Figure  33.  Sih  and  Chen  [30]  model  of  crack' in  fiber  composite 


Figure  34.  Cross  section  of  fiber  and  matrix. 


and  the  volume  of  the  fiber  VR  is 


Using  the  definition  of  the  fiber  volume  fraction,  v^  is 


V 2 

R _ rtR 

vf  m v 2 

1 T (h  + 2R) 


( 


Values  for  $(1)  and  i|r(l)  of  unity  in  Equation  (59)  would  correspond  to 
the  case  of  a crack  in  an  isotropic  plate.  Because  the  values  are  con- 
siderably less  than  unity  in  Equation  (60),  it  can  be  seen  that  the 
effect  of  the  fiber  reinforcement  is  to  lower  the  crack-tip  stress  inten- 
sities, resulting  in  a higher  fracture  strength  for  the  fiber  composite. 

The  S — theory  predicts  the  initial  direction  6 of  crack  growth 
c o 

from  Equation  (53) . Except  for  the  special  case  of  tension  normal  to  the 

crack  plane,  Qq  will  not  be  zero.  As  noted  in  Section  3.1,  the  crack 

in  a unidirectional  laminate  normally  grows  parallel  to  the  fibers, 
seemingly  at  variance  with  the  Sc~  theory.  However,  it  must  be  kept  in 

mind  that  0q  is  only  the  initial  crack  direction;  the  crack  can  turn 

after  its  initial  increment  of  extension.  For  Scotchply  1002,  Sih  and 
Chen  argue  that  the  fracture  process  is  brittle,  that  there  is  a distinct 
point  of  instability  on  the  load-deflection  curve,  and  that  most  of  the 
stored  energy  is  dissipated  at  the  point  of  incipient  fracture  where  the 
crack  extends  by  a small  amount  A a.  In  other  words,  it  is  the  initial 
increment  of  crack  extension  A a which  causes  fracture  and  the  remaining 
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energy  required  for  the  actual  separation  of  the  plate  into  two  pieces 
is  small.  Of  course,  it  is  only  the  initial  crack  direction  which  makes 
an  angle  of  0q  with  the  main  crack  surface. 

This  argument,  if  valid,  is  only  valid  for  brittle  behavior.  The 
crack  is  contained  in  a layer  of  matrix  material  whose  half  thickness  h 
is  only  0.18R.  With  no  more  room  than  this,  the  discussed  amount  of 
initial  crack  growth  Aa  at  an  angle  of  0q  would  indeed  have  to  be  small. 

Would  this  small  amount  of  extension  consume  most  of  the  fracture  energy? 
Whether  the  argument  is  strictly  correct  or  not  is  perhaps  academic  if 
the  theory  gives  good  correlation  with  experimental  data,  which  it  did 
for  Scotchply  1002  [5,  6]. 

To  apply  the  theory,  expressions  for  kj  and  k^  from  Equation  (59) 

are  substituted  into  Equation  (52)  so  that  the  energy  density  factor  S is 
known  in  terms  of  0.  Equation  (53)  is  then  used  to  find  the  initial 
crack  direction  0 . For  the  present  case,  several  different  values  of 

P,  the  load  direction,  are  being  considered  so  that  for  each  value  of  p 
there  is  a corresponding  value  of  0q.  Figure  35  shows  how  the  angle  0q 

varies  with  p.  For  p = 90°,  0q  = 0,  the  crack  direction  is  along  the 

original  crack  plane,  as  expected.  For  each  value  of  p the  value  of 
0q  can  be  inserted  into  Equation  (52) . Furthermore,  let  critical  condi- 
tions apply,  i.e.,  S -> and  a ->  ac*  Then  the  following  fracture  cri- 
terion results: 

2 T 4 3 

S = a a 0.084a..  sin  p 4-  0.049a..  sin  p cos  p .... 

c c L 11  12  (61) 

+ 0.029a22  cos4  pj 


Figure  35.  Fracture  angle  as  a function  of  crack  onele  f 301 . 
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Sih  and  Chen  determined  a value  of  S from  the  data  [5,  6].  With  the 

c 

material  constant  Sc  known,  fromEquation  (61)  the  way  in  which  critical  stress 

oc  varies  with  P can  be  determined  for  several  chosen  crack  lengths. 

Figure  36  shows  the  results  for  five  selected  half  crack  lengths: 
a ■ 0.01,  0.02,  0.05,  0.1,  and  0.51  in.  The  results  from  the  experiments 
[5,  6]  for  a =>  0.51  in.  are  shown  as  dark  squares.  Except  for  small 
values  of  (3,  the  curve  for  a a 0.51  in.  agrees  well  with  the  data.  The 
curve  for  a — 0.02  in.  agrees  well  with  a for  uncracked  specimens  [24], 

represented  by  dark  triangles.  This  indicates  that  microflaws  on  the 
order  of  0.04  in.  in  length  must  have  existed  in  the  uncracked  specimens. 
Lauraitis'  microscopic  examination  [24]  had  previously  revealed  flaws  in 
the  form  of  bubbles  and  gaps  along  the  fibers.  She  had  already  inferred 
the  flaw  length  to  be  approximately  0.035  in.  which  agrees  well  with  the 
2a  = 0.04  in.  value  indicated  in  Figure  36. 
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Figure  36.  Critical  stress  versus  crack  angle  [30]. 

The  variation  of  the  critical  stress  with  crack  length  for  several 
values  of  p is  shown  in  Figure  37.  Again,  there  is  good  agreement  between 
the  experimental  data  and  the  curves. 

The  Sc  theory  appears  to  be  a versatile  and  successful  fracture  cri- 
terion. It  will  predict  noncollinear  crack  growth  for  situations  where 
a mixture  of  k^  and  k^  exist.  The  method  deserves  further  consideration. 
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Figure  37.  Variations  of  critical  stress  with  half 
crack  length  [ 30] . 


To  verify  its  usefulness,  it  should  be  applied  to  a variety  of  unidirec- 
tional composites.  Attempts  should  be  made  to  adapt  the  theory  to 
cross-ply  laminates.  From  an  engineering  viewpoint,  the  main  disadvan- 
tage in  applying  the  model  discussed  here  is  the  difficulty  of  finding 
a stresi  analysis  for  a particular  geometry,  or  loading,  keeping  in  mind 
that  the  model  requires  the  crack  to  be  contained  in  a strip  of  isotropic 
material  surrounded  by  orthotropic  material.  Certain  approximations  may 
be  possible  to  alleviate  the  problem.  Certainly,  considerable  simplifi- 
cation resulted  [30]  for  finding  <J>(1)  and  \(r(l)  when  h is  small;  h would 
be  small  for  most  composites  of  Interest.  The  form  of  kj  and  in 

Equation  (59)  suggests  a possible  approximation.  For  a finite  plate  it 
might  be  possible  to  determine  k^  and  k^  assuming  the  crack  to  be  embedded  in  a 

homogeneous  orthotropic  plate  and  then  apply  ♦(l)  and  \|r(l)  determined 
from  an  infinite  plate  solution  as  correction  factors  to  account  for 
the  crack  lying  in  an  isotropic  strip  surrounded  by  orthotropic  material. 
This  idea  is  speculation  and  has  not  been  tested. 


69 


r 


Chapter  4.  TOUGHNESS  TESTING  OF  COMPOSITES 

Chapter  3 included  some  fracture  toughness  information;  however, 
there  the  emphasis  was  on  the  fundamental  question  of  the  applicability 
of  LEFM  to  composites,  and  toughness  testing  was  an  Incidental  issue. 

For  the  testing  of  metals,  rigorous  American  Society  for  Testing  Materials 
(ASTM)  specifications  [26]  exist.  To  obtain  valid  data,  certain  require- 
ments must  be  met  with  regard  to,  among  other  things,  specimen  thickness, 
crack  length,  shape  of  the  fracture  surface,  and  linearity  of  the  load- 
displacement  curve.  A test  that  meets  all  of  the  ASTM  standards  is  then 
referred  to  as  a valid  test  and  the  valid  fracture  toughness  is  denoted 
by  KjC*  Until  it  is  shown  that  the  ASTM  test  requirements  are  satisfied, 

the  fracture  toughness  is  referred  to  as  Kp,  the  candidate  fracture 

[toughness . Because  in  the  case  of  composites  there  are  no  standards  for 

verifying  toughness  val  idity  or  status, the  Kp  notation  has  been  generally 

adopted  to  denote  a composite's  fracture  toughness.  The  word  valid, 
when  it  is  used  in  connection  with  a composite's  fracture  toughness, 
means  only  that  the  value  is  believed  to  be  geometrically  independent, 
a material  constant;  it  does  not  mean  that  ASTM  requirements  were  met; 
this  is  impossible  by  definition. 

, 

In  this  chapter,  a sample  of  results  which  address  some  of  the 
unique  problems  of  toughness  testing  of  composites  is  included.  Emphasis 
will  be  on  general  material  behavior,  techniques,  and  methods  rather 
than  on  individual  values  of  fracture  toughness  for  any  given  material. 

4.1  Characteristics  of  Fracture  Toughness 

Fracture  toughness  testing  technology  was  developed  for  use  on 
metals.  It  is  well  to  consider  the  various  ways  in  which  the  toughness 
behavior  of  composites  differs  from  that  of  metals.  One  very  strong 
contrast  between  the  toughness  behavior  of  composites  and  metals  is  the 
difference  in  the  relationships  between  the  ultimate  strength  and  the 
fracture  toughness  for  the  two.  The  ultimate  tensile  strength  and  tough- 
ness of  a composite  are  directly  related  to  the  tensile  behavior  of  the 
filaments.  This  is  indicated  in  the  curves  of  Figure  38  for  a number  of 
laminate  configurations.  The  ultimate  strength  increases  as  the  per- 
centage of  plies  in  the  load  direction  increases.  Because  boron  is 
significantly  stronger  than  graphite,  the  boron  curve  is  steeper.  Fig- 
ure 38  (b)  shows  that  in  addition  to  the  ultimate  strength  the  fracture 
toughness  also  increases  as  the  percentage  of  load-oriented  plies  increases . 
The  fracture  toughness  and  ultimate  strength  are  directly  related.  As 
the  ultimate  strength  increases,  the  fracture  toughness  increases.  This 
is  in  marked  contrast  to  the  behavior  of  metals  which  usually  exhibit  an 
inverse  relationship  between  the  ultimate  strength  and  fracture  toughness. 
Metals  of  high  tensile  strength  are  normally  brittle. 
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PERCENTAGE  OF  PLIES  IN  THE  0°  DIRECTION 


PERCENTAGE  OF  PLIES  IN  THE  0°  DIRECTION 
(b) 
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Figure  38.  Behavior  of  strength  and  toughness  with  increasing 
load-oriented  plies  [22]. 
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In  the  design  of  a composite  fracture  test,  the  procedures,  methods, 
and  requirements  [26]  for  the  testing  of  metals  are  frequently  used  as 
guidelines.  For  example,  the  compact  tension  specimen  and  the  three- 
point  bend  specimen,  the  two  standard  specimens  for  metals  have  been  used 
considerably.  But  other  requirements  are  not  amenable.  For  example, 
fatigue  crack  cannot  normally  be  generated  in  a composite.  Various 
procedures  have  been  used  for  finding  P , the  critical  load  in  a test. 

Some  researchers  use  the  maximum  load,  some  use  the  pop-in  value,  and 
others  use  the  procedure  of  Reference  26.  Additionally,  the  mechanism 
of  crack  extension  is  far  more  complex  than  it  is  for  metals.  In  multi- 
angle laminates,  the  tendency  for  the  crack  to  follow  the  fibers  and 
the  resulting  delamination  between  plies  results  in  a very  irregular 
fracture  surface. 


In  a unidirectional  laminate,  the  tendency  for  the  crack  to  follow 
the  fibers  results  in  a combination  of  Mode  I and  Mode  II  extension 
which  is  usually  not  separable  in  the  toughness  test.  The  object  of 
toughness  testing  is  to  obtain  for  a given  material  a value  of  the  criti- 
cal stress  intensity  factor  which  will  be  independent  of  the  specimen 
dimensions.  This  raises  questions  relative  to  minimum  permissible  speci- 
men thickness,  crack  length  (both  of  which  for  metals  are  controlled  by 
the  crack-tip  plastic  zone  size),  and  specimen  width. 

4.2  Size  and  Shape  Effects 


If  the  critical  stress  intensity  factor  is  a material  constant, 
then  it  should  be  possible  to  determine  a value  which  is  independent  of 
the  specimen  geometry;  i.e.,  a number  of  different  specimen  shapes  should 
all  produce  the  same  value  of  the  fracture  toughness,  as  long  as  certain 
minimum  size  effects  with  regard  to  specimen  width,  crack  length,  etc. 
are  met.  Tests  can  be  run  using  specimens  of  varying  size  and  shape  to 
see  if  it  is  possible  to  obtain  geometry-independent  toughness  values 
and  to  determine  what  minimum  specimen  sizes  produce  anomalies  in  measured 
toughness  values. 

Cruse  and  Osias  [15]  ran  such  a set  of  experiments.  The  results 
from  three-point  bend  specimens  were  compared  with  those  from  CCT  speci- 
mens. Figure  39  shows  the  geometries.  The  thicknesses  were  considerably 
different  for  the  two  specimen  shapes:  0.07  in.  for  the  CCT  specimen 
and  0.35  in.  for  the  bend  specimen.  Moreover,  tests  were  run  for  three 
different  sizes  of  CCT  specimens,  having  scale  ratios  of  one,  two,  and 
five;  all  dimensions  except  the  thickness  were  changed.  The  material 
was  a graphite  epoxy  in  a laminate  configuration  of  (0/±45)g.  Five  beam 

specimens  and  three  center  cracked  specimens  were  tested  for  each  configu- 
ration. The  cracks  were  sharpened  ultrasonically  to  a width  of  approxi- 
mately 0.010  in.  The  critical  load  was  found  by  the  ASTM  method  [26]; 

the  critical  stress  intensity  factor,  Kq,  was  calculated  using  the 
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isotropic  formulas  [ 1] . The  results  are  shown  graphically  in  Figure  40. 
As  the  crack  length  is  increased,  K.  appears  to  approach  a constant  value 


asymtotically  (approximately  36,000  psi-in. 
the  fracture  toughness  for  the  laminate. 


) , which  might  be  called 


In  metals  the  elastic  crack-tip  stress  field  is  assumed  to  control 
fracture  if  the  crack-tip  plastic  zone  size  is  small  compared  to  the 
crack  length.  In  the  case  of  a fiber  composite,  there  is  a crack-tip 
damage  zone  which  is  analogous  to  the  plastic  zone  in  isotropic  metals. 
Cruse  and  Osias  [ 15]  suggest  that  with  the  fibers  perpendicular  to  the 
crack  the  damage  zone  size,  rQ,  can  be  estimated  from  a formula  similar 

to  that  for  metals,  namely 


/K  \2 

- _L  (-2.) 

2*  \aiu/ 


where  is  the  candidate  fracture  toughness  and  cx^  is  the  fiber  direc- 
tion ultimate  stress.  If  rQ  is  small  compared  to  the  crack  length,  then 

the  crack-tip  elastic  stress  field  would  be  expected  to  control  the 
fracture  behavior.  If  rQ  becomes  too  large,  then  the  usual  elastic  stress 

field  characterized  by  the  stress  intensity  factor  is  no  longer  an  accurate 
representation  of  the  crack-tip  stresses  and  geometrically  independent 
results  for  would  not  be  expected.  In  Figure  40,  the  ratio  of  half 

crack  size  to  damage  zone  size  is  approximately  fourteen  for  the  longest 
crack  and  approximately  six  for  the  shortest.  Perhaps  a better  way  to 
view  this  is  by  estimating  the  relative  sizes  of  the  damage  zone  and 
the  elastic  singularity.  If  it  is  assumed  that  the  elastic  singularity 
is  on  the  order  of,  say,  a/5  in  size,  then  the  sizes  of  the  elastic 
singularities  for  the  CCT  specimen  are  0.10,  0.04,  and  0.02  in.  Only 
for  the  large  specimen  (a  ■ 0.5  in.)  is  the  elastic  singularity  sub- 
stantially larger  (factor  of  three)  than  the  damage  zone.  The  large 
damage  zone  size  indicates  that  due  to  effective  crack  blunting, 

should  not  be  too  sensitive  to  crack  sharpness.  In  fact,  Broutman  [33] 
has  indicated  that  the  critical  notch  tip  radius  in  a composite  is  at 
least  an  order  of  magnitude  greater  than  that  in  most  metallic  materials. 

Another  size  effect  which  Cruse  and  Osias  [ 15]  considered  is  that 
due  to  the  fiber  ply  thickness.  In  metals  it  is  assumed  that  inhomo- 
geneities are  small  compared  to  the  size  of  the  stress  singularities; 
otherwise,  the  crack  tip  may  experience  a mixture  of  modes.  In  composites, 
this  implies  that  the  ply  thickness  should  be  small  compared  to  the  singu- 
larities. The  ratio  of  singularity  size  to  ply  thickness  (approxi- 
mately 0.0064  in.)  varied  from  more  than  15  for  the  large  CCT  specimen 
to  more  than  3 for  the  smallest. 
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CRACK  HALF-LENGTH,  a (in.) 


Figure  40.  Fracture  toughness  versus  crack  length  for  the  three- 
point  bend  and  CCI  specimens  [ 15] . 
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Results  In  Figure  40  indicate  that  there  was  no  strong  thickness 
effect.  There  is  excellent  agreement  between  the  0.35-in.  thick  bending 
specimen  and  the  0.074-in.  CCT  specimen.  If  debonding  beyond  the  crack 
tip  is  small  and  not  affected  by  thickness,  then  no  significant  thickness 
effects  are  expected  [ 15] . 

While  more  work  will  be  required  before  precise  limits  can  be  placed 
on  specimen  sizes  and  crack  lengths,  it  appears  [15]  that  there  are  no 
thickness  effects  down  to  0.070  in.,  and  that  inhomogeneous  effects  due 
to  ply  thickness  are  not  significant  for  crack  lengths,  a,  greater  than 
0.5  in.  It  is  suggested  that  the  same  conclusions  should  apply  to 
boron/epoxy  laminates. 

A number  of  papers  not  discussed  here  may  be  of  interest  to  stu- 
dents. The  CCT  specimen  has  been  used  by  a number  of  researchers  in 
addition  to  Osias  and  Cruse  [15]:  by  Zimmer  [34]  for  a bidirectional 
composite;  by  Beaumont  and  Phillips  [35]  who  observed  the  effect  of 
varied  fiber  matrix  bond  strengths;  and  by  Phillips  [36]  who  compared 
the  critical  strain  energy  release  rate  with  work-of-fracture  measure- 
ments. The  toughness  of  a carbon-carbon  composite,  a material  for  ther- 
mal protection  on  reentry  vehicles,  was  investigated  by  Guess  and  Hoover 
[37]  using  a beam  in  three-point  bending. 

4.3  Tensile  Edge-Cracked  Specimens 

The  single-edge-notched  (SEN)  and  DEC  or  notched  specimens  (Figure 
41)  were  used  by  Mandell  [21]  to  Investigate  the  toughness  of  various 
graphlte/epoxles.  In  calculating  K^,  Mandell  used  a special  K-calibra- 

tion  derived  form  a finite  element  program.  The  crack,  cut  with  a dia- 
mond saw,  had  a thickness  of  approximately  0.011  in.  The  specimen  length 
was  three  times  the  width. 

The  laminates  tested  were  separated  into  two  classes,  brittle  and 
nonbrittle.  The  load  displacement  curve  for  the  nonbrittle  laminates 
deviated  strongly  from  linearity  prior  to  catastrophic  fracture.  Follow- 
ing the  nonlinearity,  the  final  fracture  was  still  sudden  with  a sharp 
drop  in  load.  The  inelastic  behavior  was  exhibited  by  a general  cracking 
of  the  specimen  throughout  the  cross  section  with  cracks  propagating 
parallel  to  the  fibers  within  each  ply.  Laminates  which  displayed  the 
nonbrittle  behavior  had  most  of  the  plies  oriented  at  45°  or  greater  to 
the  load  direction  (Table  7) . Fracture  mechanics  is  inappropriate  for 
the  nonbrittle  laminates. 

For  the  first  two  laminates  in  Table  7,  fracture  tests  were  run 
using  the  SEN  specimen  for  varying  crack  lengths:  a ■ 0.30,  0.60  and 
0.88  in.  Figure  42  shows  the  results.  Data  scatter  prevents  firm  con- 
clusions but  no  significant  variation  of  Kq  with  crack  length  is  apparent 

for  the  three  crack  lengths  Investigated. 
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SINGLE-EDGE-NOTCHED  DOUBLE-EDGE-NOTCHED 

Figure  41.  SEN  and  DEC  tensile  fracture  specimens. 
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TABLE  7.  BRITTLE  AND  NONBRITTLE  GRAPHITE/EPOXY  PLY 
CONFIGURATIONS  [21] 

I Brittle  Behavior  \ 


(90/0/90/0) a 
(45/0/-45/90) t 

<±15>s 

(±30) 8 


Nonbrlttle  Behavior 


(±45), 


(±45), 

1 

( ±60) ( 
(±75). 


Table  8 contains  results  from  the  DEC  specimen.  The  two  laminates 
of  Figure  42  are  included.  The  Kp  values  in  Table  8 from  the  DEC  speci- 
mens are  significantly  lower  than  the  SEN  values  in  Figure  42.  These 
are  directly  compared  on  Figure  42  where  the  DEC  results  are  plotted  as 
open  circles.  Mandell  [21]  suggests  that  this  difference  may  be  related 
to  a statistical  variation  in  Kp  or  possibly  to  specimen  gripping  problems. 

It  is  interesting  to  note  that  if  the  K calibration  for  the  DEC  specimen 
is  low  as  suggested  by  Figure  16  that  the  values  of  KQ  for  that  specimen 
would  also  be  depressed.  " 

The  SEN  tension  specimen  was  also  used  by  Olster  and  Woodbury  [22] 
to  find  the  fracture  toughness  of  several  glass,  boron,  and  graphite 
epoxies.  The  plate  width  was  0.750  in.  and  the  edge  crack  length  was 
0.230  in.,  a bit  short  for  valid  data  (Section  4.2).  Results  are  shown 
in  Table  9.  Several  replica  tests  were  run  for  each  configuration.  Even 
though  the  crack  length  was  below  the  desired  minimum,  the  data  in 
Table  9 serve  as  good  indicators  of  the  effect  of  lay-up  configuration 
on  fracture  toughness.  Two  values  of  critical  k^  are  given:  one 

corresponds  to  the  onset  of  fracture  and,  while  it  was  not  specifically 
stated,  the  other  probably  corresponds  to  the  maximum  load.  It  is  felt 
that  the  onset  value  should  be  used  in  residual  strength  calculations. 

This  procedure  is  at  the  least  conservative.  Also,  the  increase  in  k^ 
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TABLE  8.  CALCULATED  VALUES  OF  FRACTURE  TOUGHNESS  ASSUMING  BRITTLE 
BEHAVIOR,  DOUBLE-EDGE-NOTCHED  SPECIMENS  (GRAPHITE  EPOXY)  [21] 


Ply 

Configuration 

Specimen 

Geometry  (in.) 

Fracture 

Toughness* 

Kq  (ksiyin.) 

Width 

Crack  Length 

(90/0/90/0) 

2.0 

0.25 

28.5 

2.0 

0.50 

27.6 

(45/0/-45/90)  a 

2.0 

0.25 

30.3 

(±15)8 

1.50 

0.25 

47.1 

(±30) s 

1.50 

0.25 

27.9 

(±45)s 

1.50 

0.25 

15.1 

(±60)g 

1.50 

0.25 

6.7 

(±75) a 

1.50 

0.25 

2.8 

♦Average  of  two  or  more  tests. 


from  the  onset  value  to  the  higher  value  signifies,  if  not  complete 
failure,  at  least  the  incurrence  of  more  damage  to  the  structure.  Thus, 
using  the  lower  value  is  consistent  with  the  damage  containment  philosophy. 

4.4  The  Compact  Tension  Specimen 


The  compact  tension  specimen  has  been  used  extensively  in  the  frac- 
ture toughness  testing  of  metals.  Slepetz  and  Carlson  [25]  have  also 
investigated  this  specimen.  Figure  43,  for  the  testing  of  composites. 
They  tested  both  unidirectional  and  cross  ply  (equal  number  of  plies  in 
orthogonal  directions)  configurations  for  a glass/epoxy  and  two  types 
of  graphite/epoxies.  Material  properties  are  given  in  Table  10.  The 
laminate  thicknesses  were:  eight  plies  of  glass  and  twelve  plies  of 
graphite.  It  was  necessary  to  constrain  the  cross-ply  specimen  with 
side  plates  to  prevent  buckling.  Because  the  compliance  technique  was 
used,  the  critical  strain  energy  release  rate,  g^,  was  determined  rather 

than  Kq.  Of  course,  the  two  are  related  by  Equation  (24)  for  colinear 

crack  growth. 


To  determine  the  suitability  of  using  analytical  rather  than 
experimental  calibration,  Slepetz  and  Carlson  [25]  conducted  compliance 
tests.  The  technique  involved  applying  a load  to  a specimen  which 
contained  an  initial  slot.  From  the  load  deflection  curve,  the  compliance 
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TABLE  9.  FRACTURE  TOUGHNESS  VALUES  FOR  SEVERAL  BORON,  GRAPHITE, 
AND  GLASS  EPOXIES  [22] 


Material 

At  Onset  of  Fracture 
(ksi-in.1/2) 

kQl 

Fracture  Toughness 
(ksi-in.  ^2) 

Boroiv'Epoxy 

(0/±45/0)s 

15.2 

21.1 

(02/±45/0/90)s 

25.8 

32.6 

(0/±602/0)s 

14.8 

18.2 

Graphite/Epoxy 

(0/±45/0)s 

21.6 

27. 52 

(02/±45/90)s 

13.5 

20. 82 

(02/±603/02)t 

17.4 

26.3 

Glass/Epoxy 

(0/±45/0)s 

39 

3 

59.2 

Notes: 

1.  Recall  Kq  m k^. 

2.  Values  are  computed  but  are  not  considered  the  true  notch 
toughness  because  mode  of  failure  was  by  interply  separa- 
tion with  lack  of  uniform  crack  path. 

3.  Values  are  computed  but  are  not  considered  the  true  notch 
toughness  because  specimens  exhibited  longitudinal  crack- 
ing and  finally  failed  as  a tension  specimen  having  a 
reduced  cross  section. 


(reciprocal  stiffness)  was  determined  for  the  given  crack  length.  The 
load  was  increased  until  crack  extension  occurred.  Then  the  crack  typically 
arrested  due  to  load  drop-off  assolcated  with  fixed  grip  conditions.  The 
new  crack  length  was  measured  and  the  specimen  was  unloaded  and  reloaded 
to  obtain  a new  compliance  for  the  new  crack  length.  This  process  was 
continued  until  the  crack  had  propagated  approximately  807.  of  the  dlsr 
tance  across  the  specimen. 
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Figure  43.  Compact  tension  specimen  configuration  [25] 
TABLE  10.  PROPERTIES  OF  LAMINATES  TESTED  [25] 


Material 

a 

E2 

(lb/in. 

G12 

2 X 106) 

v12 

Fiber 

Vo lume  (%) 

1002  S -glass/epoxy 

Unidirectional 

6.9 

HI 

0.28 

55 

Cross  Ply 

4.7 

■M 

0.14 

MOD  1-5208  graphite/epoxy 

N 

Unidirectional 

19.9 

1.11 

mSm 

0.32 

40 

Cross  Ply 

12.0 

12.0 

0.021 

MOD  11-5208  graphite/epoxy 

Unidirectional 

19.5 

0.96 

1.02 

0.32 

60 

Cross  Ply 

12.1 

12.1 

0.95 

0.035 

The  strain  energy  release  rate  can  be  determined  in  terms  of  the 
compliance,  C,  [25]: 

dC 

8 " 2B  da 

where  F is  the  load,  B is  the  thickness  and  a is  the  crack  length. 


critical  value,  gc,  is  the  value  at  which  the  crack  grows  unstable  before 
being  arrested.  This  occurs  at  a load  value,  F£,  which  was  taken  to  be 
the  maximum  load  in  the  cycle  of  crack  growth. 
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The  compliance  results  were  Interesting.  In  Figure  44  for  the 
unidirectional  S -glass/epoxy , the  experimental  compliance  (open  circles) 
is  compared  with  results  from  a finite  element  program  (dashed  lines) . 

The  experimental  compliance  is  considerably  lower  (specimen  stiffer) 
than  the  analytical  results  indicate.  The  apparent  reason  for  this  is 
that  the  unidirectional  specimens  exhibited  fiber  bridging  behind  the 
advancing  crack  tip,  holding  the  fracture  surfaces  together  somewhat 
and  increasing  the  stiffness.  To  study  this  effect,  Sleptz  and  Carlson 
[25]  machined  away  these  fibers  after  each  cycle  of  crack  growth  in  tests 
on  several  S-glass  specimens.  The  compliance  values  for  the  machined 
cracks  are  shown  in  Figure  44  as  solid  circles.  The  values  for  the 
machined  cracks  are  higher  than  the  finite  element  results. 

The  compliance  for  the  graphite  cross-ply  specimens  is  shown  in 
Figure  45.  In  contrast  to  the  unidirectional  specimens,  the  finite 
element  compliance  values  were  considerably  below  those  observed  experi- 
mentally. Slepetz  and  Carlson  [25]  suggest  that  some  splitting  and 
delamination  of  the  cross  plies  may  have  increased  the  compliance  in 
Figure  45.  In  view  of  the  wide  variance  of  the  finite  element  solution 
with  the  experimental  compliance,  the  decision  was  made  to  use  the 
experimental  compliance  in  finding  g^. 

The  fracture  toughness  values  for  the  unidirectional  S-glass/ epoxy 
are  shown  in  Figure  46  for  the  case  of  a = 90°  (fibers  parallel 
to  the  crack).  The  values  obtained  from  the  specimens  with  natural 
cracks  depend  upon  the  crack  length  up  to  approximately  30%  of 
the  specimen  width.  After  that,  g^  levels  off  at  a mean  value 

of  approximately  7.6  lb/in.  The  variation  with  length  is  believed  to 
have  been  due  to  the  fiber  bridging  action  already  discussed.  Though 
considerable  scatter  is  present,  the  values  obtained  with  the  machined 
cracks  indicate  no  dependency  upon  crack  length.  The  average  for  the 
machined  crack  was  3.18  lb/in.  The  results  for  the  natural  crack  are 
considerably  higher  than  for  the  machined  crack  indicating  that  the  fiber 
bridging  contributed  significantly  to  the  toughness. 

The  unidirectional  results  for  both  types  of  graphite/epoxy  are 
shown  in  Figure  47,  for  a =*  0 and  a ■ 90°.  The  values  for  a “ 0°  could 
be  obtained  only  by  cutting  a deep  side  groove  in  a cross-ply  specimen 
such  that  only  two  plies  remain.  These  two  plies  were  normal  to  the 
plane  of  the  crack.  The  mean  value  for  a ■ 0 was  approximately  355  lb/in. 
This  is  more  than  two  orders  of  magnitude  larger  than  for  a ■ 90°.  The 
value  is  also  considerably  higher  (factor  of  two)  than  the  value  of  117 
lb/in.  in  Table  6,  obtained  by  Konish  et  al.  [7].  The  a ■ 90°  specimens 
exhibit  a great  deal  of  scatter  but  no  significant  variation  with  crack 
length.  The  averages  of  values  for  the  Mod  I and  II  materials , respectively, 
were  2.30  lb/ in.  and  1.88  lb/ in.,  again  roughly  twice  the  value  of  0.943  lb/ in. 
found  by  Konish  et  al.  [7].  Though  the  materials  used  by  Konish,  et  al.  [7] 
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0 0.2  0.4  0.6  0.8  1.0 
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Figure  44.  Compliance  versus  crack  length  curves  for 
S -glass/epoxy  specimens  with  machined  and  natural 
cracks  [ 25] . 

and  Slepetz  and  Carlson  [25]  were  both  graphite  epoxies,  they  may 
exhibit  different  fracture  toughness  values  because  of  fiber  volume 
ratios,  etc. 

The  cross-ply  results  (plies  at  0°  and  90°)  for  graphite  are  shown 
in  Figure  47.  The  mean  values  of  g c were  124  lb/in.  and  117  lb/in.  for 

Mod  1 and  Mod  II,  respectively.  For  the  cross  plies  oriented  at 
45°,  g£  values  were  not  found  because  the  crack  propagated  along  45° 

lines  in  a zig-zag  fashion  making  the  compliance  technique  questionable. 
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0 0.2  0.4  0.6  0.8  1.0 

•/W 

Figure  45.  Compliance  versus  crack  length  curves  for  Mod  I 
and  Mod  II  cross-ply  specimens  (a  « 0 and  90°)  [25]. 

The  S-glass  cross  plies  proved  to  be  so  tough  that  no  crack  could 
be  made  to  propagate.  Instead,  a system  of  superficial  cracks  developed 
over  a broad  area  of  the  notch  tip  and  the  load  displacement  curve 
became  highly  nonlinear  in  a way  characteristic  of  gross  plastic  deforma- 
tion of  a metal. 

The  tests  on  the  unidirectional  materials  for  a ■ 30°,  45°,  60°,  and 
90°  could  not  be  analyzed  by  the  compliance  technique  because  the  crack 
propagated  along  the  fibers  under  a combined  mode  condition. 
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4.5  Toughness  of  Aluminum  Matrix  Composites 

Thus  far,  the  discussion  of  epoxy  matrix  composites  has  indicated 
that  the  critical  value  of  K is  a material  constant.  Perhaps  it  has 
been  implied  that  this  might  be  true  for  aluminum  matrix  materials  as 
well  as  for  epoxy  ones.  One  very  strong  difference  in  the  behavior  of 
the  two  materials  is  the  ductility  exhibited  by  the  aluminum.  Thus,  it 

becomes  necessary  to  inquire  to  what  extent,  in  general,  the  toughness 

will  depend  upon  such  things  as  the  matrix  yield  strength,  fiber  diameter, 
and  fiber  volume  ratio.  Some  aluminum  matrix  materials  are  included  in 
the  survey  by  Chang  [38].  Recently,  Wright  et  al.  [39,  40]  investigated 
the  R-curve  and  J -integral  approaches  to  fracture  toughness  for  boron- 
aluminum.  The  bulk  of  the  discussion  here  will  rest  on  work  carried  out 
by  Hoover  [41].  He  addressed  the  preceding  questions  by  conducting  a 
test  program  on  two  series  of  unidirectional  boron-aluminum  (B-Al)  com- 
posites. To  investigate  the  effects  of  fiber  volume  fraction  and  fiber 
diameter,  three  composites  with  an  1100  Aluminum  (Al)  matrix  were  used. 

The  parameters  of  the  three  were  20  v/o  5.6  mil  B-1100  Al,  48  v/o  5.6 

mil  B-1100  Al  and  49  v/o  8.0  mil  B-1100  Al.  The  properties  are  given 

in  Table  11. 

TABLE  11.  FOUR-POINT  BEND  TEST*  RESULTS  FOR  B-1100  AL  COMPOSITES  [41] 


20  v/o 

5.6  mil  B-1100  Al 

48  v/o 

5.6  mil  B-1100  Al 

49  v/o 

8.  ’ mil  £-1100 

Composite  Yield  Strength 
(ksi) 

13.6 

23.2 

21.9 

Standard  Deviation  (ksi) 

1.23 

0.53 

1.86 

Composite  Ultimate  Strength 
(ksi) 

88.4 

221.3 

199.7 

Standard  Deviation  (ksi) 

7.9 

18.8 

11.8 

Elastic  Modulus  (psl) 

16.15  X 106 

30.4  X 106 

28.7  X 106 

Strain-to-Fracture  (X) 

0.875 

0.90 

0.84 

Standard  Deviation  (X) 

0.05 

0.04 

0.01 

*A11  testa  were  conducted  on  specimens  0.400  In.  wide  and  0.203  in.  thick  with  a 
minor  span  of  3.0  In.  and  a major  span  of  7.5  In.  The  strain  was  measured  with 
a strain  gage  on  the  tensile  surface.  A minimum  of  three  samples  of  each  mate- 
rial were  tested. 
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A second  series  of  tests  was  run  on  specimens  taken  from  a plate  of 
48  v/o  B-6061  A1  and  subjected  to  three  different  heat  treatments.  The 
resulting  specimen  had  three  different  matrix  and  composite  yield 
strengths.  Properties  are  given  in  Table  12.  These  were  used  then 


TABLE  12.  FOUR- POINT  BEND  TEST*  RESULTS  FOR  B-6061  AL  COMPOSITES  [41] 


48  v/o 

48  v/o 

48  v/o 

B-6061-01 

B-6061-F2 

B-6061-T63 

Composite  Yield  Strength 

(ksi) 

17.1 

37.2 

99.0 

Standard  Deviation  (ksi) 

0.52 

4.5 

2.9 

Composite  Ultimate  Strength 
(ksi) 

228 

235 

244.3 

Standard  Deviation  (ksi) 

10.4 

4.6 

10.0 

Strain-to-Fracture  (%) 

0.94 

0.98 

0.97 

Standard  Deviation  (%) 

0.04 

0.01 

0.01 

*A11  tests  were  conducted  on  specimens  0.400  in.  wide  and  0.208  in. 
thick  with  a minor  span  of  3.0  in.  and  a major  span  of  7.5  in. 

The  strain  was  measured  with  a strain  gage  on  the  tensile  surface. 
A minimum  of  three  samples  of  each  condition  were  tested. 


1. 

2. 

3. 


Annealed  for  2 hours  at  420°  C and  furnace  cooled. 

As  fabricated  condition. 

Solution  treated  for  30  min.  at  490°C,  ice-water  quenched 
and  aged  for  18  hours  at  155 °C. 


to  investigate  the  effect  of  matrix  yeild  strength  on  toughness.  Hoover 
[ 41]  used  the  rule  of  mixtures  to  calculate  the  effective  yield  strengths 
for  the  three  materials.  Matrix  yield  strengths  were:  5.7  ksi  for  the 
0 condition,  12.4  ksi  for  the  F condition,  and  33.0  ksi  for  the  T6  condi- 
tion. 


The  specimens  were  loaded  in  four-point  bending.  The  specimens  were 
0.400  in.  deep,  0.200  in.  thick  with  a span  of  1.57  in.  The  starter  notch 
varied  from  0.120  to  0.330  in.  The  notch  root  radius  varied  from  0.0005 
to  0.005  in.  The  fibers  were  normal  to  the  notch. 
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Figure  48  shows  a load-deflection  curve  for  the  5.6  mil  B-1100  Al. 
The  curves  exhibited  substantial  nonlinearity  and  a pop-in  load  value, 

Pc,  characterized  by  a small  dip  in  the  curve.  The  pop-in  load  P was 

used  to  calculate  K^,  using  the  formula  for  isotropic  materials.  It 

should  be  noted  that  the  pop-in  load  is  considerably  higher  than  the 
critical  load  which  would  be  found  from  the  graphical  method  of  ASTM. 
Using  the  pop-in  value  for  all  tests,  however,  should  lead  to  a nt 

consistent  comparison  for  evaluating  the  effects  of  matrix  ductility, 3 
etc.  Table  13  shows  the  values  found  by  Hoover  for  various  crack 

lengths  in  the  B-1100  Al  composites.  It  can  be  &r-»n  that 


0 0.02  0.04  0.06  0.08  0.10 

MIO-SPAN  DEFLECTION  (in.) 

Figure  48.  Typical  load-deflection  behavior  of  a 
5.6  mil  B-1100  Al  composite  [41]. 

Kq  has  a constant  value  for  all  crack  lengths  except  for  a 0.32  in. 

For  that  case,  the  width  of  the  uncracked  portion  of  the  beam  was  only 
0.08  in.,  i.e.,  a free  boundary  was  very  close  to  the  crack  tip.  The 
values  of  found  for  the  smallest  crack  length,  a « 0.12  in.,  were 

in  agreement  with  the  others.  That  such  small  crack  lengths  would  yield 

a constant  value  of  K is  surprising.  For  metals,  ASTM  requires  that  the 

" 2 
crack  length  and  specimen  thickness  both  exceed  2.5  where 

is  the  critical  stress  intensity  factor  and  a is  the  yield  strength. 

In  the  present  case,  If  this  criterion  is  evaluated  using  the  composite 
yield  strength  from  Table  11  then  a crack  length  and  specimen  thickness 
of  several  inches  results.  This  criterion  appears  unrealistic  for  B-Al 
composites,  since  Hoover  found  constant  values  for  all  crack  lengths 

except  the  longest.  Recent  tests  by  Ads it  and  Waszczak  [42]  indicate 
also  that  for  B-Al  is  a constant.  Their  results,  based  primarily  on 

two  specimens  (a  tensile  specimen  with  a center  hole  and  a DEC  specimen) , 
indicate  a minimum  required  flaw  size  considerably  higher  than  Hoover's. 
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TABLE  13.  TOUGHNESS  DATA  FOR  B-llOO  AL  COMPOSITES  [41] 


Sample 


20  v/o  5 . 6 ml  1 B- 1 100  A1 


Crack  Length  Crack  Length  to 
a,  (In.)  Width  Ratio,  a/w 


Fracture  Toughness 

KQ(k8iytn.) 


48  v/o  5.6  mil  B-1100  A1 


49  v/o  8.0  mil  B-1100  A1 


0.279 

0.315 


0.118 

0.156 

0.198 

0.238 

0.281 

0.318 


0.115 

0.158 

0.195 

0.242 

0.276 

0.324 


0.498 

0.618 

0.698 

0.788 


0.295 

0.390 

0.495 

0.595 

0.703 

0.795 


0.288 

0.395 

0.488 

0.605 

0.690 

0.810 


Their  curve  shown  inFigure  49  monotonically  approaches  a constant  value  of 
Kq  at  a flaw  size  of  approximately  0.3  in.  as  contrasted  with  Hoover's 

minimum  crack  size  of  approximately  0.12  in.  Their  curve  contains  the 
combined  influence  of  two  specimens —one  with  a hole  and  one  with  DECi 
From  the  viewpoint  of  fracture  toughness  testing,  a hole  does  not  simu 
late  a crack.  Thus,  the  hole  specimen  may  have  unduly  influenced  the 
variation  of  with  flaw  size.  The  main  point  is  that  KL  for 
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Figure  49.  of  unidirectional  boron/aluminum  as  a function 
of  flaw  size  [42]. 


unidirectional  B-Al  does  appear  to  be  a material  constant  independent  of 
crack  size  and  specimen  dimensions. 

Table  13  shows  the  influence  of  the  fiber  volume  fraction.  The 

1/2 

toughness  for  the  48  v/o  (approximately  28  ksi/in.  ' ) is  approximately 

1/2 

607o  greater  than  the  20  v/o  (approximately  17  ksi-in.  ' ) . The  effect 
of  fiber  diameter  is  also  shown  in  Table  13.  There  is  no  apparent 
difference  in  the  toughness  of  the>5.6-  and  8.0-mil  materials. 

The  effect  of  crack  tip  plasticity  as  influenced  by  heat  treatment 
of  the  B 6061  A1  is  shown  in  the  load-deflection  curves  of  Figure  50. 

Both  the  annealed  (0  condition)  material  and  the  as-fabricated  (F  condi- 
tion) material  exhibit  considerable  nonlinearity  prior  to  the  pop-in 
load,  as  is  usually  characteristic  of  metals  with  substantial  crack  tip 
plasticity.  The  T6  material  has  a higher  pop-in  load  and  a more  nearly 
linear  load-displacement  curve.  The  T6  heat  treatment  increases  the 
yield  strength,  thus  decreasing  the  extent  of  crack  tip  plasticity  so 
that  more  nearly  linear  behavior  results.  The  heat  treatment  did  not 
appear  to  change  the  value  of  deflection  at  which  pop-in  occurred. 

The  resulting  fracture  toughness  values  for  the  three  heat  treatments  is  shown 
in  Table  14. 

Hoover  points  out  that  when  the  load  is  applied  to  the  composite, 
work  is  required  at  the  crack  tip  to  stretch  the  fibers  elastically  and 
to  deform  the  A1  matrix  elastically  and  plastically.  The  crack-tip  work 
is  thus  partitioned  between  the  fibers  and  matrix.  The  yield  strength 
of  the  matrix  influences  the  partitioning  between  the  fibers  and  matrix. 
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Figure  50.  The  effect  of  heat  treatment  on  the  load-deflection 
response  of  B-6061  A1  composite  [41]. 

Because  the  matrix  yield  strength  did  not  affect  the  deflection  at  which 
crack  initiation  occurred,  the  difference  in  toughness  is  due  to  the 
difference  in  energy  partitioning  at  the  crack  tip.  The  material  with 
the  highest  yield  strength  (T6  treatment)  absorbs  more  energy  to  the 
matrix,  resulting  in  a higher  toughness. 


Hancock  and  Swanson  [43]  have  reported  a fracture  toughness  value 
1/2 

of  33.7  ksi-in.  ' for  30  v/o,  4 mil,  B-6061 -T6,  Al.  As  already  noted  for  a 
higher  volume  fraction,  the  toughness  should  be  greater.  The  value  in 
Table  15  for  a volume  fraction  of  48  is,  in  fact,  larger  by  approximately 
30%.  Adsit  and  Waszczak  [42]  reported  a toughness  value  for  5.6  mil, 
B-6061,  Alin  F condition;  the  fiber  volume  fraction  was  not  reported. 

1/2 

The  toughness  reported  by  them  of  106.7  ksi-in.  seems  too  high  in 

comparison  with  the  value  of  approximately  34  ksi-in.1^2  in  Table  14 
for  the  F condition  material. 


TABLE  14.  TOUGHNESS  DATA  FOR  48  V/O  5.6  MIL  B-6061  AL  COMPOSITES  [41] 


Sample 

Crack  Length 
a,  (in.) 

Crack  Length  to 
Width  Ratio,  a/w 

Fracture  Toughness 

Kc(kstyin.) 

B-6061-  0-1 

IB2I 

0.273 

25.6 

0-2 

BaSH 

0.353 

27.5 

0-3 

MEM 

0.478 

27.6 

0-4 

B8SH 

0.593 

28.5 

0-5 

UULLJU 

0.678 

26.0 

0-6 

0.313 

0.783 

18.9 

B-6061-  F-l 

0.120 

0.300 

33.5 

F-2 

0.158 

0.395 

33.9 

F-3 

0.197 

0.493 

36.7 

F-4 

0.242 

0.605 

34.3 

F-5 

0.265 

0.663"' 

30.3 

F-6 

0.325 

0.813 

21.1 

B-6061-  T6-1 

0.116 

0.290 

45.0 

T6-2 

0.150 

0.375 

45.1 

T6-3 

0.192 

0.480 

43.6 

T6-4 

0.226 

0.565 

44.2 

T6-5 

0.293 

0.733 

44.2 

T6-6 

0.311 

0.778 

25.9 
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Chapter  5.  FRACTURE  OF  COMPOSITES  CONTAINING  STRESS  CONCENTRATIONS 

The  application  of  fracture  mechanics  to  a structural  component 
of  a given  material  requires  a knowledge  of  that  material's  fracture 
toughness,  whether  obtained  from  the  literature  or  from  a fracture  tough- 
ness testing  program  as  discussed  in  Chapter  4.  This  chapter  also 
discusses  fracture  toughness,  but  in  light  of  three  very  specific  appli- 
cations relating  to  composites  containing  various  stress  concentrations. 

Certain  fiber  composites  with  notches  display  peculiar  fatigue 
behavior:  surviving  specimens  of  a fatigue  test  may  show  a residual 

strength  equal  to  or  greater  than  a specimen  which  has  not  been  fatigue 
loaded.  Moreover,  the  tensile  static  strength  of  laminated  composites 
containing  circular  holes  depends  upon  the  hole  size.  This  is  not  pre- 
dicted by  the  classical  stress  concentration  factor.  These  two  features 
will  be  considered  in  this  chapter. 

5.1  Inherent  Flaw  Model 

The  behavior  of  fiber  composites  with  stress  raisers  presents  a 
marked  contrast  with  that  of  metalB.  Typically,  fatigue  loading  of 
notched  specimens  fails  to  produce  fatigue  cracks.  The  residual  strength 
of  the  notched  survivors  of  a fatigue  test  is  almost  the  same  as  for 
unnotched  specimens.  Moreover,  the  fracture  of  such  specimens,  in  the 
absence  of  a distinct  fatigue  crack,  is  nonetheless  brittle.  Furthermore, 
static  tests  of  specimens  with  circular  holes  have  shown  that  the  laminate 
strength  varies  with  hole  size,  even  though  the  stress  concentration 
factor  does  not.  Waddoups  et  al.  [44]  presented  the  data  in  Table  15, 
which  show  the  residual  strength  of  notched,  graphite/epoxy  laminates 
after  fatigue  loading.  The  residual  strength  of  the  notched  laminates 
is  nearly  the  same  as  the  unnotched  laminate.  The  residual  strength  of 
the  notched,  fatigued  laminates  exceeded  the  static  strength  of  those 
not  fatigued. 

Tables  16  and  17  present  data  by  Waddoups  et  al.  [44]  based  on 
tests  performed  at  General  Dynamics.  These  tables  show  the  effect 
of  hole  size  on  laminate  strength.  As  hole  size  increases,  the  static 
strength  decreases. 

Waddoups  [44]  introduced  a model  to  explain  the  preceding  behavior. 
The  model  consists  of  a hole  from  which  two  intense  energy  regions 
radiate  (Figure  51) . The  Intense  energy  regions  are  modeled  as  cracks 
which  have  a small  but  finite  length  a.  The  length  a is  taken  to  be  the 
characteristic  length  of  an  inherent  flaw.  If  it  is  now  assumed  that 
the  isotropic  stress  intensity  factor  applies  to  the  geometry  of  Fig- 
ure 51,  then  the  Bowie  solution  [1]  or  [3]  for  cracks  emanating  from 
a hole  can  be  used.  The  Bowie  solution  is  written  as 
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TABLE  15.  [ 0/90]  GRAPHITE/EPOXY  FATIGUE*  [44] 


Specimen** 
Diameter  (in.) 

Average  Static  Strength 
(psi) 

Average  Residual  Strength 
After  5 X 10^  cycles  (psi) 

No  stress  concentration 

83,500 

77,200 

0.063  Hole 

68,900 

77,700 

0.063  Hole  + 

0.018  X 0.004  Notch 

72,200 

74,300 

0.063  Hole  + 

0.078  X 0.004  Notch 

58,300 

79,400 

*This  study  was  conducted  under  an  in-house  experimental  program  at 
General  Dynamics  Convair  Aerospace  Division  Fort  Worth  Operation  by 
E.  L.  McKague  and  R.  J.  Stout. 

**Coupon  specimens  1.0  in.  wide  and  9.0  in.  long. 


TABLE  16.  SMALL  HOLE  DATA  SUMMARY  [44] 


Specimen**  Hole 

Diameter  (in.) 

Static  Strength  (psi) 

Control 

67,240 

0.062 

45,000 

0.031 

51,500 

0.015 

60,900 

*Thi8  work  was  sponsored  by  the  Air  Force 
Materials  Laboratory  under  Contract  F33615- 
69-C-1494. 

**Coupon  specimens  were  1.0  in.  wide  and 
9.0  in.  long. 


where  f (a/R)  is  a correction  factor  accounting  for  the  hole  boundary. 

Let  critical  conditions  apply,  -»  k , and  d ->  cr  , where  k is  the 

x y c y 

critical  stress  intensity  factor  and  cr^  is  the  critical  stress; 
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TABLE  17.  LARGE  HOLE  DATA  SUMMARY*  [44] 


Specimen**  Hole 
Diameter  (in.) 

Static  Strength  (psi) 

Actual 

Corrected  for  Finite  Width 

Control 

76,000 

76,000 

1.0 

26,600 

27,900 

2.5 

15,900 

22,800 

3.0 

13,250 

23,000 

,fThis  work  was  sponsored  by  the  Air  Force  Materials 
Laboratory  under  Contract  F33615-69-C-1494. 

**Coupon  specimens  were  5.0  in.  wide  and  38.0  in.  long. 


I 


Because  a is  assumed  to  be  small,  the  free  edge  correction  factor  is  unity. 
Assuming  that  is  a material  constant,  the  control  specimen  failed  at 


a critical  stress  of  a so  that 

o 


ao'/*® 


From  Equations  (64)  and  (66) 


<*c  \R 


thus,  the  ratio  of  the  failure  stress  of  a control  specimen  and  the 

failure  stress  crc  of  a ‘notched  specimen  is  given  by  the  Bowie  correction 

factor  if  a,  the  inherent  flaw  length,  is  known.  Figure  52  shows  the 
ratio  0Q/ac  in  terms  of  the  hole  radius  R for  four  assumed  values  of  a. 


HOLE  RADIUS  R (In.) 

Figure  52.  Parametric  study  of  critical  stress  and  flaw  length  [44]. 

The  failure  stresses  fromlables  16  and  17  are  also  plotted  there.  It  can  be 
seen  that  a value  of  a * 0.04  in.  fits  both  sets  of  data  points  with  reasonable 
accuracy. 
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Based  on  the  strength  of  a control  specimen  and  a notched  specimen 
with  the  Waddoups  model  it  is  possible  to  predict  the  strength  of  other 
notched  laminates  of  similar  geometry.  The  method  is  a two-parameter 
one;  the  two  parameters  are  and  a.  Equation  (67)  is  used  to  find  the 

function  f(a/R)*  With  f(^/R)  known,  the  Bowie  solution  is  used  to  find 
the  ratio  a/R.  For  the  notched  test  specimen,  R is  known  and  with  a/R 
known,  the  inherent  flaw  length  a can  be  calculated.  With  a known,  then 
either  Equation  (64)  or  Equation  (66)  can  be  used  to  find  K^,  the  critical 

stress  intensity  factor.  With  and  a known.  Equation  (64)  can  be  used 

to  predict  the  critical  stress  a of  other  notched  laminates.  A numerical 

c 

example  will  clarify  the  procedure.  If  the  data  for  hole  radii  less  than 
0.5  in.  in  Table  16  are  considered  and  the  strengths  of  the  first  two  specimens 
are  substituted  specimens  into  Equation  (67) , the  result  is 


f(*\  m 12  . 67x240 

\RJ  a 45,000  L' 


495 


From  Bowie's  results  [1]  or  [3] 


t * °'89 


so  that. 


a - 0.89R  - 0.89^|^)  - 0.0276  in. 

Kq  can  now  be  determined  from  Equation  (64): 

Kq  - ocy«a  f(|j  - 45,000 ,/jr( .0276)  1.495 


K 


Q 


19,800  psi-in. 


1/2 


Knowing  K^  and  a,  the  inherent  flaw  length,  the  static  strength  of  the 
other  two  notched  specimens  can  be  determined.  From  Equation  (64) 


The  specimen  with  the  0.031  in.  diameter  hole  is  considered: 


R- 


0.0155 
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a _ 0.0276  . 7. 

R 0.0155  A‘,B 

and  from  the  Bowie  solution. 


(Note:  Some  of  the  numbers  here  differ  slightly  from  those  in  Reference  44.) 
From  Equation  (68): 

crc  54,200  psi 


This  agrees  reasonably  well  with  the  measured  value  of  51,500  psi.  For 
the  specimen  with  a hole  of  0.015  in.,  ac  is  found  to  be  60,200  psi  as 

compared  to  the  measured  value  of  60,900  psi. 

The  data  in  Table  17  are  now  considered  for  hole  diameters  of 
1.0  to  3.0  in.  Using  a for  the  control  specimen  and  a for  the  1-in. 
hole,  ° c 

f(a\  . fo  „ 76,000  _ 
ryR  J oc  27,900 

from  the  Bowie  solution 


f - 0.100 


so  that 

a - 0.100  - 0.0500  in. 

and  from  Equation  (64) 


M ■ *■ 


Kq  - 30,100  psi-in. 
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For  the  specimen  with  the  2.50-in.  hole, 

a.  2^. „.M0 


and 


f “3-13 


then  o ■ 24,300  psi 
c 
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which  compares  well  with  the  measured  value  of  22,800  psi.  Similarly, 

for  the  laminate  with  a 3.0-in.  hole,  a is  found  to  be  23,900  psi  as 

c 

compared  to  a measured  value  of  23,000  psi. 

The  values  found  for  the  inherent  flaw  lengths  of  0.0276  and  0.0500  in. 
bracket  the  nominal  value  of  0.0400  in.  noted  before  from  Figure  52. 

The  physical  basis  for  an  Inherent  flaw,  a,  has  been  enforced  some- 
what by  Cruse  [ 27] . Cruse  modeled  circular  notches  with  equivalent  ideal 
cracks.  The  equivalent  crack  length  was  obtained  by  calculating  stress 
intensity  factors  for  the  notches  based  on  the  similarity,  over  a given 
region  of  the  transverse  stresses  ahead  of  a notch  and  a crack.  The 
equivalent  half-crack  length  was  greater  than  the  hole  radius  by  an 
amount  which  was  close  to  the  inherent  flaw  length.  The  inherent  flaw 
length  calculated  from  Equation  (63) , was  nearly  the  same  for  the  two 
laminates  (0^/±45)  and  (0/±45) . This  implied  the  possibility  that  the 

inherent  flaw  size  is  real  and  relatable  to  the  microstructure. 

The  Waddoups  model  has  been  discussed  by  Brinson  and  Yeow  [45]  who 
compared  the  model  stresses  with  measured  failure  stresses  on  several 
graphite  epoxy  laminates  containing  holes.  The  comparison  was  not  direct; 
it  was  necessary  to  adjust  the  data  for  the  finite-width  plates  before 
comparison  with  the  theory  which  is  for  infinitely  wide  plates.  The  cor- 
relation between  theory  and  experiment  was  better  for  unidirectional 
laminates  than  for  angle  ply  laminates. 

It  has  been  pointed  out  by  Kanninen  et  al.  [46]  from  a philosophical 
viewpoint  that  the  Waddoups  model  is  essentially  empirical,  a two  param- 
eter empirical  correlation  of  limited  test  data.  Nevertheless,  the 
model  may  afford  a useful  method  of  predicting  the  strength  of  notched 
structural  components.  It  needs  further  verification  for  a number  of 
laminate  materials  and  configurations. 

5.2  Stress  Fracture  Criteria  for  Notched  Laminates 

The  fracture  of  laminates  containing  holes  has  been  investigated  by 
Whitney  and  Nuismer  [47,  48]  using  two  stress  criteria.  While  it  is  not 
a fracture  mechanics  approach,  their  method  assumes  the  existence  of  an 
inherent  flaw.  By  applying  their  criteria  to  the  case  of  a crack,  they 
derive  a basic  relationship  for  the  effect  of  crack  length  on  fracture 
toughness. 


A hole  of  radius  R in  an  orthotropic  laminate  subjected  to  a remote 
stress  of  a parallel  to  the  y-axis  is  considered.  The  normal  stress 
o along  the  x-axls  in  front  of  the  hole  is 


(69) 


(x,  o)  - | (2  + 


$ [<f)6-7(f)8]) 


x > R 


where  K^.  is  the  orthotropic  stress  concentration  factor  for  an  infinite 
plate  and  is  determined  from 


" 1 +/c22  KC11C22 


c , C11C22  - C12  \ 

12  2 C22  ) ' 


The  constants  are  the  in-plane  laminate  stiffness  coefficients  of 

Equation  (12).  In  terms  of  the  effective  elastic  moduli  for  the  laminate 
which  are  the  same  form  as  Equations  (3) , but  indicated  by  an  overline 
here.  Equation  (70)  becomes, 

h ■ 1 +A{.fA-  vn)  + fr:  <71> 


For  the  isotropic  case  K^,  ■ 3,  Equation  (69)  reduces  to  the  first  three 
terms  which  are  the  exact  isotropic  distribution. 
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For  the  sake  of  simplicity,  the  remaining  discussion  will  be  based  on  this 
equation  for  an  isotropic  or  quasi-isotropic  laminate  rather  than  the  more 
complex  expression  in  Equation  (69).  If  Equation  (72)  is  plotted  in  the 
usual  manner  as  a function  of  the  nondimens ional  variable,  x/R,  then  the 
results  obviously  are  independent  of  the  hole  size.  In  terms  of  physical 
distance  from  the  edge  of  the  hole  (x-R) , however,  the  distribution  will 
be  different  for  different  hole  sizes.  For  example.  Figure  53  shows  the 
distribution  for  two  cases:  R - 1.0  and  R ■ 0.1.  It  seems  reasonable 
that  the  strength  of  the  laminate  with  the  larger  hole  would  be  lower. 

For  the  larger  hole,  a greater  volume  of  material  is  subjected  to  high 
stress,  increasing  the  probability  that  an  inherent  flaw  will  exist  in 
a region  of  high  stress,  thus  lowering  the  strength  of  the  plate.  From 
Figure  53  it  also  appears  that  the  smaller  hole  has  a better  opportunity 
to  redistribute  the  stress  resulting  in  a higher  strength  for  the  plate 
with  the  small  hole.  These  two  observations  suggested  to  Whitney  and 
Nuismer  that  a criterion  based  on  the  stress  magnitude  at  a point  — in  this 
Instance,  the  edge  of  the  hole— was  not  reasonable  and  that  the  stress 
distribution  over  some  characteristic  length  must  be  considered. 


Figure  53.  Stress  distribution  for  a hole  in  an  infinite  isotropic 
plate  [ 47] . 


They  chose  two  models,  each  involving  a characteristic  length 

denoted  by  d and  a . The  first  model  assumes  that  failure  occurs  when 
o o 

the  normal  stress  at  a fixed  distance,  dQ,  from  the  edge  of  the  hole 

reaches  a value  of  oq,  the  unnotched  tensile  strength  of  the  laminate. 

The  distance  d^  is  assumed  to  be  a material  constant;  it  represents 

the  distance  over  which  the  material  must  be  highly  stressed  to  find 
a flaw  sufficiently  large  to  initiate  fracture.  The  criterion  in 
equation  form  is 


a (o,  x) 

yy 


r + d 


Using  Equation  (72)  results  in  the  ratio  of  the  notched  to  unnotch  strength 


(*♦«!♦  < 


where  oN  is  the  notched  strength  and  |^  ■ R/(R  + dQ) . For  small  holes. 


0 and 


1 and  Equation  (74)  predicts  no  strength  reduction; 


whereas  for  large  holes,  | 


and  cr^/c 


1/3,  so  that  the  strength 


reduction  is  governed  by  the  classical  stress  concentration  factor.  Thus 
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Equation  (74)  gives  the  correct  limits.  The  appropriate  value  of  d0 

must  be  determined  from  test  data.  Whitney  and  Nuismer  [47]  compared 
Equation  (74)  using  varying  values  of  dQ  to  test  data  on  a quasi-isotropic 

glass/epoxy  of  different  hole  radii.  Their  results  are  shown  in  Figure  54. 
The  values  of  dQ  are  the  same  as  those  found  by  Waddoups  [44]  for 

boron/epoxy  using  the  Bowie  crack  model:  0.03  in.  fits  the  small  hole 
data,  0.05  in.  fits  the  large  hole  data,  and  0.04  in.  is  a good  compromise 
value  to  represent  all  data  with  reasonable  accuracy. 


Figure  54.  Stress  reduction  in  quasi-isotropic  glass/epoxy 
laminates  due  to  the  presence  of  a circular  hole,  point 
stress  criterion  [47]. 

A second  approach  developed  by  Whitney  and  Nuismer  [47]  requires 
calculating  the  average  stress  over  some  distance  aQ,  also  taken  to  be 

a material  constant  Independent  of  laminate  construction  and  stress 
distribution.  Physically,  this  distance  is  assumed  to  represent  a length 
over  which  failure  has  occurred  with  subsequent  stress  redistribution. 
Failure  is  assumed  to  occur  when  the  average  stress  reaches  the  unnotched 
tensile  strength,  oq.  In  equation  form, 

, R+a 

r 1°  ayy  (x’  0)d*  " °o  • <75> 

Substituting  Equation  (72)  into  Equation  (75)  and  doing  the  integration 
results  in 
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2«  - 52) 

— * 7 — 24\  (76> 

° (2  - S2  - *2) 

where  = R/(R  + aQ).  Equation  (76)  is  compared  with  glass/epoxy  data 
in  Figure  55.  The  values  of  aQ  shown  are  0.1,  0.15,  and  0.2  in.,  with 
the  value  of  0.15  in.  giving  the  best  comparison  with  data. 


tt/4  SCOTCHPLY 
0o  « 45.3  ksi 


•c- 0.15  in. 


•Q  - 0.2  in. 


ay  v “ °o 


•o-0.1  in. 


(X-R) 


Figure  55.  Stress  reduction  in  quasi-isotropic  glass/epoxy  laminates 
due  to  the  presence  of  holes,  average  stress  criteria  [47]. 

The  data  in  Figures  54  and  55  are  too  limited  to  conclude  whether 
the  point  stress  or  the  average  stress  criteria  is  most  accurate.  Nuismer 
and  Whitney  [47]  have  compared  the  two  models  for  a number  of  laminate 
configurations  of  both  glass/epoxies  and  graphite/epoxies.  The  compari- 
son, based  on  the  general  anisotropic  results  in  Equation  (69), was  incon- 
clusive due  to  data  scatter.  Brinson  and  Yeow  [45]  compared  the  two 
models  with  test  data  on  several  graphite  epoxy  laminates.  After  adjust- 
ing the  finite-plate  data  to  the  infinite  plate  case,  the  comparison 
was  rather  inconclusive.  As  a design  tool,  either  of  the  two  models 
seems  reasonably  accurate. 

The  two  models  of  Whitney  and  Nuismer  [47]  also  provide  insight 
into  the  question  of  dependency  on  crack  length.  A straight  crack 

of  length  2a,  in  an  infinite  anisotropic  plate  loaded  remotely  by  a 
stress  normal  to  the  crack  surface  is  considered.  If  the  origin  of  the 
x-y  coordinate  system  is  located  at  the  crack  center,  the  normal  stress 
a ahead  of  the  crack  tip  is  approximated  by 


(77) 


yy  sriT(X  - a) 

This  is  the  usual  fracture  mechanics  stress.  For  this  case,  ■ ajiia.. 

Equation  (77)  is  accurate  at  points  within,  say,  a/10  of  the  crack  tip. 
Obviously  for  short  cracks,  small  a,  the  region  of  applicability  is  very 
narrow.  Therefore,  Whitney  and  Nuismer  proposed  using  the  exact  expression 


f ita  (x  - a ) 


/ 2 2 
'x.  - a 


Equations  (77)  and  (78)  are  compared  in  Figure  56  for  a » 0.1  and  1.0  in. 
The  stress  in  Equation  (77)  as  a function  of  (x  - a)  , is  independent  of 
the  crack  length  and  gives  the  same  curve  for  both  crack  lengths;  this  is  shown 


I 1, 


0.04  0.06 
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Figure  56.  Stress  distribution  for  a crack  in  an  infinite 
anisotropic  plate. 

as  a dashed  line.  Equation  (78)  shows  a considerable  difference  in  the 
stress  for  the  two  crack  lengths.  For  the  long  crack,  the  approximate 
and  exact  expressions, Equations  (77)  and  (78),  show  good  agreement.  It 
is  apparent  that  if  failure  is  assoicated  with  a fixed  characteristic 
length  such  as  dQ  or  aQ,  constant  for  all  crack  lengths;  then  for  short 

cracks.  Equation  (78)  should  be  used. 


Now  the  exact  stress  Equation  (78)  is  used  in  both  failure  criteria 
(changing  R to  a)  in  Equations  (73)  and  (75) . For  the  point  stress  and 
average  stress  criteria,  respectively,  results  are 


(79) 


where 


a + d 


and 


fa  A23 
V i + 5, 


(80) 


where 


a + a 


Equations  (79)  and  (80)  give  the  strength  reduction  for  a crack  according 
to  each  failure  criteria.  Both  indicate  no  strength  reduction  for  short 
cracks  and  vanishing  strength  as  a becomes  large. 

To  observe  the  crack  length  size  effect  on  fracture  toughness,  K^, 
it  is  noted  that  = a^^/ica  and  Equations  (79)  and  (80)  are  rewritten  as 

Kq  = aQ/« a(l  - ?*)  (81) 


for  the  point  stress  criteria  and 

1 


K = a 
Q o 


Hl  - *4: 
J (i  + %,.) 


(82) 


for  the  average  stress  criteria.  For  long  crack  lengths.  Equations  (81) 
and  (82)  asymptotically  approach  the  respective  constant  values  of 


K, 


a J2nd 
o v o 
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(83) 


- 

1 


and 


K = cr 

Q o 


/F 


(84) 


which  result  from  the  approximate  stress  given  by  Equation  (77). 

Figures  57  and  58  show  the  crack  length  effect  for  both  criteria  com- 
pared with  graphite/epoxy  data.  The  figures  show  that  no  constant  value 
of  Kq  should  be  expected  below  a crack  length,  a,  of  approximately  0.4  in. 

This  compares  well  with  the  recommended  minimum  crack  size  noted  in 
Chapter  4 based  on  Figure  40. 


Figure  57 . Critical  stress  intensity  factor  as  a function  of  crack 
length  for  graphit^/epoxy  laminates,  point  stress  criterion  [ 47] . 

5.3  Fracture  Due  to  Ballistic  Impact 

The  hazard  of  ballistic  impact  is  one  which  must  be  considered  in 
the  design  of  military  aircraft.  Two  quantities  are  of  immediate  inter- 
est: the  threshold  strength  and  the  residual  strength.  The  threshold 
strength  is  the  applied  load  at  which  the  structure  will  immediately 
fail  upon  impact  by  a given  type  of  projectile  impact.  The  residual 
strength  refers  to  the  maximum  load  the  structure  will  still  sustain 
subsequent  to  projectile  damage.  Analytically,  the  determination  of 
threshold  strength  is  a dynamic  problem  influenced  by  time  dependent 
material  response  and  the  nature  of  the  target-projectile  interaction. 
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Figure  58.  Critical  stress  intensity  factor  as  a function  of  crack 
length  for  graphitu/epoxy  laminates,  average  stress  criterion  [47], 


On  the  other  hand,  the  determination  of  residual  strength  of  a damaged 
component  is  normally  a static  problem,  involving  only  the  fracture 
toughness  of  the  material  and  a characterization  of  the  projectile 
damage . 

For  a number  of  lay-ups  of  both  graphite/epoxy  and  bororv/epoxy, 

' Olster  and  Woodbury  [ 22]  experimentally  determined  the  threshold  strength  and 

residual  strength  of  panels  subjected  to  30-caliber  (0.30-in.  diameter) 
armor  piercing  projectile  impact.  Figure  59  schematically  shows  the  loaded 
specimens  in  the  experimental  setup.  The  measured  residual  strengths  of 
graphite/epoxy  panels  for  various  lay-ups  varied  from  approximately 
61%  to  737.  of  the  true  ultimate  tensile  strength.  The  threshold  strengths 
were  slightly  lower,  varying  from  approximately  51%  to  65%  of  the  tensile 
strength.  For  boron/epoxy,  the  residual  strength  varied  from  approxi- 
mately 527.  to  657.  of  the  tensile  strength.  Again,  the  threshold  strength 
was  only  slightly  less,  approximately  92%  of  the  residual  strength.  The 
measured  residual  strengths  for  the  various  lay-ups  are  shown  in  Table  18. 

Because  the  threshold  strength  appears  to  be  only  slightly  less  than 
the  residual  strength,  predictions  of  the  residual  strength  also  provide 
a rough  approximation  to  the  threshold  strength.  The  advantage  of  this 
observation  is  that  residual  strength  predictions  are  much  easier  to 
make  than  are  threshold  predictions.  To  make  residual  strength  calcula- 
tions, the  character  of  the  expected  damage  and  the  fracture  toughness 
of  the  material  must  be  known. 
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Figure  59.  Fracture  mechanics  model  of  a laminate  panel  damaged  by 
ballistic  impact  of  a 30-caliber  armor  piercing  projectile. 


TABLE  18.  COMPARISON  OF  THE  PREDICTED  AND  MEASURED  RESIDUAL 
STRENGTH  OF  A BALLI STI CALLY  DAMAGED  PANEL 


Laminate 

Configuration 

Fracture 

Toughness 

kc 

(ksiyin.) 

Residual  Strength 

Predicted 

Measured 

1 

Measured 

(ksi) 

Percent  of 
True  Ultimate 
Strength 

UW 

Boron  Epoxy 

(0/±45/0)s 

21.1 

46.8 

44.8 

52.6 

(0z/±45/0/90)g 

32.6 

72.6 

74.9 

65.6 

(0/±602/0)s 

18.2 

40.5 

41.2 

60.6 

Graphite  Epoxy 

(0/±45/0)b 

21.6* 

48.0 

51.9 

73.0 

(02/±45/90)g 

13.5* 

30.0 

36.0 

63.1 

(02/± 603/02)T 

17.4* 

38.6 

32.9 

61.4 

*At  onset  of  failure. 
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Olster  and  Woodbury  [22],  in  their  experiments,  characterized  the 
projectile  damage  as  a hole  from  which  cracks  radiate.  The  hole  diameter 
is  approximately  the  same  as  the  bullet  diameter.  The  hole-edge  cracks 
were  measured  (although  somewhat  crudely)  to  be  approximately  0.04  in. 
in  length.  They  used  the  same  crack  length  for  both  boron  and  graphite. 
Now  this  problem  (Figure  59)  is  the  same  as  the  Waddoups  [44]  model. 

In  fact,  it  is  interesting  to  note  that  the  crack  length  is  the  same 
as  the  inherent  flaw  length  assumed  by  Waddoups. 

Although  some  of  the  values  are  suspect,  the  fracture  toughness  for 
the  subject  laminates  determined  by  Olster  and  Woodbury  [22]  were  already 
discussed  and  are  shown  in  Table  9;  the  appropriate  values  are  repeated 
here  in  Table  18.  The  stress  intensity  factor  solution  needed  for  the 
problem  is 

kj  - Yaa1//2  (85) 

where  the  values  of  Y are  given  by  Figure  22  for  the  various  laminates. 
However,  for  a crack  length  of  0.04  in.,  the  isotropic  solution  gives  a 
good  approximation  to  all  of  the  laminates.  Therefore,  Y « 2.24  is  used 
so  that 

kj  » 2.24  aa1>/2  (86) 

It  should  be  noted  that  the  definition  of  k^  being  used  here  is 

= Kj/^/n.  In  Equation  (86)  let  k^  assume  its  critical  value  for  each 

respective  laminate  and  a “ 0.04  in.;  then  o,  the  residual  strength, 
is  easily  calculated. 

These  calculated  values  of  residual  strength  are  compared  with  the 
measured  residual  strength  in  Table  18.  There  is  good  agreement  for 
boron/epoxy.  The  agreement  for  graphite/epoxy  is  fair.  It  must  be 
explained  that  different  definitions  of  fracture  toughness  were  used  for 
the  two  materials.  Two  sets  of  values  were  given  in  Table  9,  one  [22] 
corresponds  to  the  fracture  toughness;  the  other  corresponds  to 
the  value  of  k at  the  onset  of  fracture.  These  two  definitions 
were  used  for  boron  and  graphite,  respectively.  This  is  _an 
unfortunate  ambiguity  caused  by  the  lack  o.t  a working  definition  of 
exactly  what  constitutes  fracture  (i.e.,  a definition  in  terms  of  load 
displacement  behavior  such  as  that  for  metals)  during  a toughness  test. 

To  contrast  the  behavior  of  the  two  materials  under  discussion  here,  it 
was  noted  [ 22]  that  during  the  toughness  tests  the  boron  materials  frac- 
tured in  a manner  characteristic  of  a valid  test  with  very  little  de lami- 
nation whereas  the  graphite  materials  suffered  considerable  crack-tip 
delamination  before  final  fracture.  Although  the  vigor  of  this  example 
is  compromised  by  the  ambiguous  use  of  the  fracture  toughness  values,  it 
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nevertheless  serves  to  illustrate  methods  of  strength  calculations  and 
the  potential  usefulness  of  making  residual  strength  predictions. 
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